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FOREWORD 


The classical differential equations for the problem of three 
bodies remain valid only if there are no collisions or other dis- 
continuities for real values of time. The equations of motion are 
not analytic when two or three of the bodies occupy coincident 
positions. In order to investigate collisions, the equations of 
motion must be made analytic by a suitable transformation of the 
independent variable. Once this transformation is carried out, the 
equations of motion are regularized. This paper is an exposition 
of Sundman’s treatment of regularization of the three body prob- 
lem. Although Sundman’s work is the basis for this paper, related 
papers and discussions have been included. To my knowledge, this 
paper is the first complete exposition of Sundman’s historic paper 
to be done in English and in vectorial form. The paper provides 
all developments in detail and leaves very little to be taken for 
granted. 


EXPOSITION OF SUNDMAN’S 


REGULARIZATION OF THE THREE BODY PROBLEM 


by 

Donald K. Yeomans 


ABSTRACT 


A complete exposition of Sundman’s regularization of the three 
body problem is given. The equations of motion and the integrals 
of motion are derived. Double real collision is investigated and 
the vector joining the center of mass of the two colliding bodies 
and the non-participating third body is found to be bounded. The 
velocity and acceleration approach infinity as the distance between 
the two colliding bodies approaches zero. The unit vector ap- 
proaches a limit near real collision. A new independent variable 
“u” is introduced which is seen to remove the singularity in the 
equation of motion for double collision. The mutual distances 
between the three bodies, along with the original independent 
variable are expanded into a power series in “u M . A lower limit 
for the strip of convergence of these solution series is determined. 
Another independent variable “w” is introduced and it is seen to 
remove all singularities for any number of double real collisions 
between any of the three bodies. A lower limit is also determined 
for the strip of convergence of the power series solution with 
respect to this variable. The convergence of the power series 
solutions is investigated and found to be extremely slow. 
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NO i FiLiVitO. 


EXPOSITION OF SUNDMAN'S REGULARIZATION OF THE THREE BODY PROBLEM 


INTRODUCTION 

The motion of the three body system is considered regular if the coordinates of the system are 
analytic functions of the independent variable. These equations of motion for the three body prob- 
lem remain analytic as long as the mutual distances between the three bodies remain greater than 
zero. A singularity is encountered when two bodies collide. The series solutions for the coordi- 
nates in the three body problem, converge only in so far as there are no singularities. The main 
body of this paper will be concerned with the removal of the singularity brought about by collision 
of two bodies. This paper gives the exposition of Karl Sundman's "Memoire sur Le Probleme Des 
Trois Corps,” although additional material and explanations are included. Sundman showed that 
the singularity of the differential equations which corresponds to a collision of two of the bodies is 
not of an essential character, and may be removed altogether by making a suitable change of the 
independent variable. The new independent variable is chosen in such a way that the differential 
equations of motion are regular and a real prolongation of the motion after collision is possible. 
The coordinates can then be specified for all values of time, whether collisions take place or not 
and a positive lower bound can be assigned to the two greater of the mutual distances. The coordi- 
nates of the three bodies, and the time are analytic functions of the new independent variable T V" 
and they can be expanded as convergent series in powers of "r" for all real values of the time. 


THE EQUATIONS OF MOTION 

If the masses of the three bodies are designated m 0 , m 1 and m 2 and their position vectors from 
an arbitrary point ”0” given by T 0 , and r 2 , the equations of motion are written: 



where A 01 , \ £ refer to distances between the masses. Adding Equations 1 we have; 


d 2 ?„ d 2 7, d 2 r 2 

0 dt 2 + m ‘ dt 2 + m 


2 dt 2 


0 


( 2 ) 



Integrating we have; 


dr 0 dr, 


dr. 


“dt" + m i dt + m 2 dt 


(3) 


A is a constant vector which provides three integrals of motion. Integrating once again we have, 


m o 7 0 + m i ^1 + m 2 r 2 “ At + B 


(4) 


B is a constant vector. We now define R as the position vector of the center of mass of the system, 
where 


_ m Q r 0 + m i F 1 +m 2 r 2 

^ ~ m Q + m 1 + m 2 


(5) 


From Equation 3 we have 


dR 

dt 


A 

M 


m Q + nij + m 2 


( 6 ) 


and from Equation 4 we have 


R 


At + B 
M 


(7) 


From Equations 6 and 7 we see that the center of mass of the system moves in a straight line with 
constant velocity (i.e., inertial system). Switching to the center of mass system we have 


R = 0 


dR 

dt 


0 


A = B = 0 

m 0 + mj r l +m 2 r 2 


0 


(8) 


If we define 


u 


f m 0 rtij f m Q m 2 f rrij m 2 


02 


(9) 


2 


We can rewrite Equations 1 as 


"o-^T = grad^U 


»l " erad ?i U 


- grad 7 U 

1 o 


dot multiplying the first Equation in 10 by dr 0 , the second by dr, and the third by dr 2 we have, 


2 d K V o + m i v, 2 +tn 2 V 2 2 ] 


d?; • d7 o + d?; • d7 i + d7^ • d7 2 = du 


finally we have by integration the energy integral 


2 L m 0 V 0 Tm l V 1 m 2 V 2 


= U + a 


where a = constant. Cross multiplying the first Equation in 1 by r 0 , the second by 7 1 and the 
third by f 2 we write, 


f m j n^Q 


f m q m j 1 3 


r o x m o dt 


A3 r o * r i + A 3 

01 02 


1 x m l dt 


f m j m q f nij m 2 r ^ x r 2 

A 3 r l * r 0 + A 3 

■ 01 12 


r 2 X m 2 dt 


f m 2 rn q f m 2 nn j r 2 ^ r j 

•V2 ?2 x 7 ° + 


3 



adding 


a 

<1 

o 

dvi 

^2 

dt + m i r l 

x dt 

+ m 2 r 2 + dt 


(15) 


Since 7 x dv/dt = d/dt (r x v) we can write 15 as 


m o^ (r o x v 0 ) +m 2 ^(r 2 xv 2 ) = 0 


(16) 


Integrating we get 


"o ( 7 0 X 7 o) + m i ( 7 1 * Vi) + m 2 ( 7 2 X V 2 ) = C 


(17) 


Rewriting Equations 8, 13 and 17 the ten integrals of motion are 


2 ( m 0 v 0 2 +m i v, 2 + m 2 v 2 2 ) = U + a 


n o ( 7 o x v 0 ) + (r, *Vj) + i" 2 (? 2 * v 2 ) - c 


m o 7 o + m, r, + m 2 r 2 - 0 


m o v o + m, v, + m 2 v 2 - 0 


(18) 


The first equation is the integral of energy, the second equation (3 integrals) is the angular mo- 
mentum integral and the last two equations (6 integrals) are referred to as the center of mass 
integrals. 

It is evident that our problem consists of nine equations in 14 of second order and hence 18 degrees 
of freedom. With the corresponding 10 integrals of motion, there remains 8 integrals for a solu- 
tion to the problem. These integrals are not known. However a reduction of the problem can be 
accomplished by using the last two equations in 18. This simplification reduces the degrees of 
freedom to 12. This reduction is accredited to Jacobi and is described below. 


JACOBI'S REDUCTION 

In the diagram (Figure 2) K is the center of mass of m 0 and m, while s represents the center of 
mass of the three body system. Our goal is to derive the equations of motion in terms of the two 
vectors r andp, thus reducing the system to 12 degrees of freedom. 


4 



recombining terms 


7 o ( m 0 + m i + m 2 ) + tHj r + m 2 (\r+^) - 0 


or 


r (m 1 + Xsn 2 ) + ~jom 2 


where 


M - nij + m 2 + m 3 


/ m 1 m 2 \ 

: \ m i + n^~nrj + ^2 


m, ( m 0 + m l + m 2 ) 

0 +m l . 


+ pi n 2 


( 21 ) 


finally 


r\ 


2 _ 

M P 


( 22 ) 


Once again we write 


Substituting Equation 22 for r 0 


Xr ■ u 


m 0 +m l 


1U 2 

IT P 




2 

M 


m 2 p 


6 


lotting 


m 0 * nij P 


r i P T ’ If P 


(23) 


(24) 


Rewriting Equation 20 


r 2 - Kr + p + r 0 


Using Equation 22 we have 


\r + p - Kr - -w 


M P 


% * m j 


TT~ p 


m 0 + m i 


M P 


(25) 


rewriting Equation 25 


p 


M _ 

m 0 + nij r 2 


(26) 


If we define 


M 

ro 2 ( m o +tn l) 


(27) 


we have from 26 


p 


Rm 2 r 2 


(28) 


7 


Collecting formulas and definitions, remember that; 


7 0 = 

m ^ 

- ^ ■ a p 

7 r 7 1 “ 7 o 

k 

II 

3 

* 

3 

o 

= 

m 2 

p* - H p 

P - e m 2 r 2 

P 

m o 

m i + m o 

? 2 = 

m 0 + nij 

M P 


g 

M 

m 2 ( m 0 + m7] 

J X 1 


At this point Sundman changes notation. The distances between the bodies are denoted 


^02 I r 2 r 0 | 


I P + K T \ 


I P ~ Mr I 


( 28 ') 


where the distance between m, andm 2 is denoted with the missing index (r 2 ). It should be emphasized 
that r 2 is not the absolute value of the vector r 2 . The notation r 2 and r are used interchangeably 
depending upon the circumstance. 


Writing the equations of motion (1) using the new notation we have: 


We write for convenience, 


d 2 ? 0 

f m, ( 


7 o) 

f m 2 | 

^ r 2 r o) 

dt 2 


r/ 


+ 

r i 3 

d 2 F, 

f m 0 

( 7 o‘ 

7 l) 

f m 2 

{ 7 2~ 7 l) 

dt 2 




+ 

r 3 
r o 

d J ? 2 

f m 0 

( 7 o~ 

? 2 ) 

f m j 

{ 7 2 " 7 l) 

dt 2 


r, 3 



r o 3 




m 2 m 0 

+ m i . 

_ 

" ^0 

kr 

+ ~M p + — 

M P - 

P 


m o 

+ m i _ 

m 2 


■ ? 1 


M P P T 

+ IT P = 

P 

~ 7 0 

- r 


r 2 = x 2 

+ i 


( 29 ) 


( 30 ) 

( 31 ) 

( 32 ) 


8 


Substracting the first Equation in 29 from the second and using 32 we write 


d 2 r 
"dt 2 



f rrij r 



+ f m 2 



( 7 2 _7 0 


r 


3 

1 


) 



(P +kr) 

r i 3 . 


So that 


d 2 r 
"dt 7 


f ( m 0 +m i) + 




(33) 


By Equation 26 


- - M _ 

P m Q +m l r 2 


differentiating we have 


d 2 p M * 2 

dt 2 m 0 + m i dt 2 


Substituting in the third equation in 29 and using 30, 31 


dt 2 


d 2 : 


m 0 +m l dt 2 m o +m l 


f m o (“ P “ ^ ) f m 1 (- p + [xr ) 
+ 


d 2 p 

It 2 


f M 


f M 


+ Hi p 1 f*L ) 

L" r . 3 r o 3 J 
* tj ) * ■ i?) 


( 34 ) 
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Rewriting Equations 33, 34 we have 



Looking at Equations 35 we have completed Jacobi's reduction. We now have 6 Equations of the 
second order for 12 degrees of freedom. Since we used the 6 integrals of the center of mass for 
the above reduction, we are left with 4 known integrals. 

We note that in a collision of masses m, and m 0 , r goes to zero while r t and r 0 remain re- 
latively large and in fact equal to each other. This can be seen by referring to Figure 3. 



With these approximations the second Equation in 35 can be written 


d 2 p _ - f M p 

dt 7 = 7 


(36) 


With this approximation the problem is reduced to the two body problem. 


AREA INTEGRAL 

The area integral will now be derived. From Equation 18 we rewrite 

dr o ^ dr, ^ dr 2 

m o r o * dt + m i r i * dt + m 2 r 2 x dt ~ c 


10 


Using Equations 22, 24, and 25 


(_ H!i -A A ( ^ -V / „ m 2 d ( _ m 2 A / m o + m i A d / m o m i A 

m o ^ Xr U P)^dt \- Kt ~ M Pj +n 'i\^--U p)^-dt[^~lA P) x di[ M Pj = c 


, 2 - dr % m 2 _ dr % m 0 in 2 _ dp m 2 2 _ dp , _ df m l m 2 _ df 

m o A - rx dt + ' v M ^ X dt + A - M rx dt + M 2 m o px dt +m i P r x dt M"~ rt° x dt 


m l m 2 _ d p m i m 2 _ dp ( m o +m i) 2 , dp _ _ 

M - x dt + ^ * dt + " ^ m 2 ' df = c 


collecting terms 


- df r , _ ap 

r x df l m o ^ + m i ^ J + P x df 


- dp 


m 2 2 m 0 nij m 2 2 (iriQ+nij) 2 


M 2 M 2 


M 2 


2 


_ df 
+ P x dt 


A-iii q m 2 fjxn j m 2 


!~] _ dp r^O m 2 m 2^| 

J + r x dt L - M ~M J 


From Equations 20, 23 the last two terms are zero and 


- dr / _ dp 

r x dt \ m o K + m i P ) + P x dt 


2 2 
m ^ ni q nil m 2 


M 2 M 2 


( m o + m iV 
M 2 


(37) 


now using 20 and 23 again 


m 0 \ 2 + m. p 2 


m 0 mj 2 + mj m 0 2 


n Q nij ( mj + m 0 ) m o m i 


( m l +m o) 2 ( m l +m o) : 


1 u, 0 


(38) 


and 


( m o +m i) 2 


M 2 M 2 M 2 


m 2 2 m Q + m l m 2 2 4 m 2 m 0 2 2m 0 m j m 2 + m 2 t 

M 2 


(™1 +m 2 + m 3)( m O +m i) 


M 2 


m 2 (m 0 + m j J 

M 


(39) 
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If we define 


h 


m o +m i 


(40) 


and from 27 


M 

m 2 (">0 + m l ) 


Using this notation, Equation 37 may be written 


_ dr _ dp _ _* 

gr x ^ + ViyO x at ” S hc 


Area integral 


(41) 


ENERGY INTEGRAL 

To derive the energy integral we start with Equation 13 


( d7 o) 2 

(dr A 2 

( d M 2 

m o \dt / +m i 

V dt ) 

+ m 2 \dt 1 


U * a 


(42) 


where now from Equation 9 


u = f i— ; r c 


n 0 m 2 \ 
“■T" / 


(43) 


Substituting Equations 22, 24, 25 into 42 we write 


/ dr m 2 dp\ / 

n o x dt M dt ) + m i V 

squaring the appropriate terms we have 


dr 
^ dt 


dp) 2 /"o^'ldpf 
dt ) + m 2\ M dt ) 


U + a 


f 

’ M 

r \ 2 2Km 2 d? 

dp m 2 / C 

p) 2 ' 


r vd 

It j + M dt 

dt + M \d 

it / J 



, /dr \ 2 

m 2 2 

(dp\ 2 

m i 

P [dlj 

M dt dt + M 2 

\dt / 


+ 


( m o +m i ) 2 

M 2 



U + a 


12 


Collecting terms 


1 

2 



k 2 + nij y 2 } 


X / dp\ 2 f-o m 2 2 
2 \ dt 1 \ M 2 


m i m 2 2 m 2 ( m 0 + m i ) 2 
+ + 

M 2 M 2 


1 dr* dp / 2 m 0 ^ m 2 2m j pm 2 ^ 

1 2 dt dt \ M ' M ) ~~ 


from Equations 20, 23, the last term on the left is zero and using Equations 38, 39, 40 and 27 we 
have 


1 

2 B 




ghU + a' 


if we define 


hgU = V 


(44) 


we have 


g 



= 2V - K 


Energy integral 


(45) 


The next few pages will be devoted to deriving some equations which will be useful in in- 
vestigating the motion of the system. Differentiating Equation 45 


dr 
2g dt 


d 2 r dp d 2 p 

+ 2h dt • 


dV dr 2 ^ dV ^ r j. dV dr o 

2 dr 2 dt + 2 drj dt 4 ^ dr Q ‘ dt 


this can be represented vectorially as 


r + v 


,*] 


d 2 r 

A , d 2 p J 

> 

) 

> 

dV ~ dV . dV\ 

8 

t+h ^ P J 

- (v 2 i +Vj j + V 0 k ) ■ 

V <^2 + J dr l +U d 'o/ 


or 


V 


d 2 r d 2 p 

a + h 

dt 2 dt 2 _ 


v * gradV 


13 



we have 


grad ? V 


h dt* - = erad?v 


The Euler Equation for a homogeneous function of the order n 


i*. * 


f = f (*i> 


Since V is a homogeneous function of order (-1) 


r • d? + P • fp = C 1 ) V 


from Equations 46 we have 


gr ' ^ + ' 1* 


adding Equations 45 and 47 we have 


which is equivalent to, 


diilill , h r«* 2 (p!) 


- d -(g? 2 +h^) = 2(V-K) 


14 


If we define 


R 2 = gr 2 + h p 2 


(50) 


then from Equation 49 


d 2 R 2 

= 2(V-K) 


(51) 


Now we assert that 


2 _ 2 r 2 


2 ^ r > r 2 

R 2 = -IT- + — + — 


(52) 


where r 0 - A J2 , r , - A Q2 , r 2 r A 0I Equation 52 is justified if it can be reduced to Equation 50; using 
Equation 28' we have 


R 2 


R 2 


_ [p 

-pr| 2 


| p + \r | 2 _ 

|r I 2 




m o 



X 

m 2 



= £i 

- 2 pp • 

r 

+ P 2 r 2 

+ eL 

+ 2 Kp 

’ 7 + \ 2 r 2 r 

2 


m o 



m j m 

2 


' 1 

1 

1 + 2 

[V + 


1 1 - - 1 

"2^ 

- p 2 

_ m o + 


j ♦ r 2 

— + 

m 0 

— + 

m i 

^J + p ‘ r l 



"o J 


by Equations 20, 23 the last bracket is zero and 


[m, +%' 

, _ 2 

m o +m i 1 

L m i m o 


_( m 0 + m l) 2 m2 _ 


or 


R 2 - P 3 


'm, + %" 

X 2 

m 2 + m 0 + m i 

m l m 0 _ 

+ r z 

( m 0 +m l) m 2 


15 



and using Equations 27, 40 we have R 2 gr 2 +h p 2 which is Equation 50 and assertion 52 is justified. 
From Equation 51 


d 

r dR] 

d 2 R 

/dR V 

dt 

I^HtJ 

= 2R dt + 

2 (dtJ 


2(V-K) 


and 


R 


d 2 R 
dt 2 



V - K 


(53) 


We now differentiate R 2 (Equation 50) to get 



- dr _ dp 

g r * dt + h p ' at 


(54) 


Squaring Equation 54 



g z r 


(at) + h2 ^ 2 (a?) + 2 h z* 


dr 

dt 


dp 

dt 


(55) 


We now define 


P = 


gh / dp dr \ 2 

l r dt - P dtj 




Since (ax b f- a 2 b 2 - (a • t>) 2 we can write 


P 



dp dr\ 2 g T 2 (dr^i 2 l - . drV 

dt -P dt ) + 77 L r ^tJ ' V dt ) 



(56) 


(57) 


Since r dr/dt z r • dr/dt and p dp/dt = p • dp/dt the last two brackets in 57 are identically zero 
and 


P 


gh / dp_ dr \ 2 
^2 \ r dt P dt j 


(58) 
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using Equations 55 and 58 we write 



from Equation 50 and Equation 45 



2V - K 


(59) 


It is evident from Equation 53 and Equation 59 that 


P 


d 2 R 

R + v 
dt 2 


or 


R 


d 2 R 
dt 2 


P - V 


(60) 


adding Equations 60 to Equation 53 we have 


2R 


d 2 R 
dt 2 




= P - K 


In an effort to express P (Equation 56) in a different form we introduce 


hR 2 
gr 2 p 2 


/- Ap gh -\ 2 

(pxTt ~v? 2c ) 


+ 


g 2 h 2 

R 2 


c 


2 


(61) 


(62) 
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from Equation 41 we denote 


ghc - c, + c 2 


_ _ _ dr 

Ci - grx^ 


where 


- = dp 

c 2 h p x ^ 


(63) 


Substituting Equations 63 into expression 62 


hR 2 
gr 2 P 


2 




“12 



J 



E 2 


squaring 


hR 2 
gr 2 p 2 



2£ - 
hR 2 Cl 


2 ^ 2 . £ 

hR 2 ° 2 R 4 


+ c 2 + 2c . 



Cj 2 + c 2 2 + 2c j 
R 2 


2 


rearranging terms 


hR 2 2hR 2 /Q 2 + hR 2 yO 4 + _1_" 

gr 2 p 2 h 2 gr 2 /J 2 hR 2 gr 2 p 2 R 4 R 2 _ 


_ I" hR 2 2 P 2 hR 2 2 p 4 + 2 1 (64) 

+ C ‘ ' ° a L gr 2 p 2 hR 2 + gr 2 p 2 R 4 + R 2 J 
the coefficient of c x • c 2 by Equation 50 is written 


2 + 2hp 2 + 2h/0 2 + 2gr 2 2_ + 2R 2 _ 

gr 2 + gr 2 R 2 + R 2 gr 2 gr 2 R 2 gr 2 gr 2 R 2 


hR 2 

Lgr 2 P : 


p_ 

R 4 


R 2 


+ c„ 


18 


therefore 64 is written 


r h p 2 

1 " 

, R2 2 hp* l] 

Lgr 2 R 2 

' R 2 . 

2 [_S r 2 /° 2 h gr 2 gr 2 R 2 R 2 J 


Since by Equation 50 


h p 2 | 1 = gr 2 + hp 2 _ 1 

gr 2 R 2 R 2 gr 2 R 2 gr 


we write 


c 


2 



R 2 


[gr 2 p 2 h 



or 



R 2 


Lgr 2 p 2 h 




by Equation 50 again we have 


c i 2 + 2 gr 2 + hp 2 -h p 2 

gr 2 ° 2 L gr 2 P 2 h 



(62') 


using Equations 63 



now this equation is equal to the expression 62* and upon further inspection it is also the last two 
terms in the original definition of P (Equation 56). Therefore we can rewrite Equation 56 as 


P 


gh 

R 2 


[ dp dr\ 2 

\ r dt ~ P dt / 


hR 2 /_ dp gh p 2 

gr V X 31 ' ~R 2 " 



g 2 h 2 

R 2 


(65) 
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Using the derived equations we shall now investigate double and triple collisions. Some important 
theorems will result. 


DOUBLE COLLISION - p REMAINS BOUNDED 
Remembering Equation 51 


d 2 R 2 
dt 2 


2(V-K) 


we shall first investigate a collision of m 0 and m, (double collision). In this case r - 0 as t -tj where 
t, is the time of collision now since V = hgU (Equation 44) and U is proportional to r' 1 , V necessarily 
goes to infinity. We write if r - 0 => V - , t - 1 , and from Equation 51 


d 2 R 2 

-^T > 0 t 1 -S 0 <t<t 1 (66) 

where tj - S 0 is an interval about tp from 66 it is evident that dR 2 /dt is increasing and we have two 
possibilities 


a) 


dR 2 

dt 


> 0 


b) 


dR^ 

dt 


< 0 


(67) 


Case a) implies R 2 is smoothly increasing with time and case b) implies R 2 is smoothly decreasing 
with time. 




Case a Case b 


Figure 4 

Either 67 a) or 67 b) could be true in a double collision. 
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We will now show that as r - o, p -p t where p, is finite. This proof will be done using the 
method of reductio ad absurdum. For a double collision, r -o and p is finite. We shall assume that 
P does not approach a limit. Since R 2 - gr 2 +hp 2 (Equation 50) it is evident that R 2 does not approach 
a limit. Since R 2 is continuous with continuous derivatives, and is without limit by assumption, R 2 
must oscillate. However dR 2 /dt, by 67 a and b, has a constant sign and R 2 cannot oscillate. We 
have a contradiction and the original assumption is incorrect. Therefore p has a limit which we 
denote p,. This completes our proof and we write Theorem 1 : 


lim p - p 
r-0 1 


From 28 T we write 


\p~P*\ 
| p + \ r | 

Irl 

as t - 1 , , r - o and 



r o = P i 


Pi 


r 2 - 0 


This trivial conclusion for a double collision was obtained from the signs of the first and 
second derivatives of R 2 . 


FOR A TRIPLE COLLISION, THE AREA INTEGRAL (C) EQUALS ZERO 
In the case of a triple collision we have only to consider case b of 67 


dR 2 

dt 


< 0 


and 


R 2 - 0 
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writing Equation 61 


d 2 R 

2R —7 + 

dt 2 



P - K 


and Equation 65 


P 


gh / dp dr \ 2 

^ l r at - P Ht ) 


hR 2 / _ d p 
gr 2 p 2 \ px a* 


ghp 2 
R 2 C ) 


R 2 


If we let the first two terms of P = F, we have 


g 2 h 2 c 2 
R 2 


+ F 


(68) 


where F >0. Substituting into Equation 61 


d 2 R /dR \ 2 g 2 h 2 c 2 

2R — + U) = 


dt 2 


+ F - K 


or 


2R 


d 2 R 
dt 2 


/dR \ 2 g 2 h 2 c 2 
+ (dt ) - g ~^~ + K = F 


(69) 


Multiplying both sides of Equation 69 by dR/dt 


dR d 2 R 
^ dt • ’ 


dt 2 


(nr - 


g 2 h 2 1 
R 2 


dR 

dt 


dR 

K dt 


dR 

F HF 


(70) 


For convenience we define 


H 



+ KR + 


g 2 h 2 c 2 


(71) 
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differentiating, we have 


dH 

dt 


dR /dR\2 
dt \dt) 


dR 
2R dt 


d 2 R 
dt 2 


K dt 


l 2 h 2 

R 2 


dR 

dt 


Since dH/dt equals the left hand side of Equation 70 we write 

dH dR dR 

dt * F ^ or dH - F af • dt 

We now stop to introduce some convenient notation; when t = t' , R = R' andH = 
t = t ' , R = R " and H = H" 


t' < t < t 


Integrating 72 and using the above notation we have 



we now consider two subcases a) and b) 

, dR 

a ) dF > 0 f° r t K t < 1 


therefore R is increasing and R' <R";from 74 


dR 

F dF dt > 0 


and H ' < H 


b) 


dR 

dt < 


for 


< t < t 


R is decreasing and R' > R” , H' > H" . Since we are considering a triple collision, 
creasing and we only consider case b) that is; 


dR 


dF < 0 


(72) 

' likewise when 

(73) 

(74) 


( 75 ) 

must be de- 
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and R^O as t - 1 x , H" < H ' , R " < R ' . Rewriting Equation 71 


H 



+ KR + 


g 2 h 2 c 2 

R 


Since R is positive 


H > KR 


g 2 h 2 c 2 

R 


and 


KR 


+ 


g 2 h 2 c 2 

W 1 


< H" 


Since H"<H' 


KR' 


g 2 h 2 c 2 

W 1 


< H' 


now as t" -tj , R" ^0 if we are to have a triple collision. If R" is then zero we have °°£H\ Since 
H* is finite, we must conclude that c = 0 . This is one of Sundman T s important theorems. That is, 
Theorem 2: For a triple collision, the area integral (c) equals zero. This is a necessary, but not 
sufficient condition. 

BEHAVIOR OF VELOCITY VECTOR NEAR DOUBLE COLLISION 

We now return to the case of the double collision in an effort to investigate the velocity near 
collision, i.e., to investigate dr/dt ast-t r For double collision r-oast-t,. In an interval 
1 1 - < t < 1 1 there exists a number b such that 


r 


o 


> b 


r 


l 


> b 


when 


r 


2 


r - 0 


(76) 
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Rewriting Equations 35, 


dt 2 


q 3 r 3j -"’ 2 ^ ri 3 + rQ 3 


~ = ~ fK+mJ -j- 3 + fm 2 p(^ - - f- — 

d 2 j3 


dt 2 


= fM 


-P[f? + + ^iTi - 71 


-/ i i 


The second of these equations can be rewritten 


d 2 jo 
dt 2 


_ (jo +\r) _ fMA. (p-fir) 

r i 2 ^ ri r o 2 r ° 


By the triangle inequality 


d 2 p\ 

fM 

p + \r 

fMX. 

p ~ pr 

dt 2 


r i 


r o 


by Equation 30 and 31 this reduces to 


and from 76 


d 2 p 

< fM 

JL + A. 

dt 2 

r 2 r 2 

J l r oJ 


d 2 p\ 

fM 


dt 2 | 

- b 2 


(77) 


in the interval 1 1 - b 0 < t < t , . 

We note that this acceleration is bounded. If we define t ' = t,-8 0 <t<t 1 then we note that 


d 1 

dt 
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where the prime indicates dip/ dt has been evaluated at time t ' . By the triangle inequality 



Since 


t 


i 


= t l ' *1 + S o 


- S 


o 


We see that I (dp/dt )' I is finite because it is evaluated before the time of collision and \ dip/ dt I 
is seen to be bounded. At a later time, I dp/dt | will be shown to have a definite limit at collision. 

Rewriting the energy integral (Equation 45) 



as t - 1 , V has been shown to approach infinity and the above discussion leaves dp/dt finite. There 
fore we realize that dr/dt , or the velocity of the colliding masses along the vector joining them, 
approaches infinity as t -t.. Thus far, we have seen that p has a definite limit, dp/dt is bounded 
and dr/dt approaches infinity as t - t t in the case of a double collision. The following section will 
be devoted to further investigation of dr/dt as the masses m 0 andmj approach collision. 

We start with r 2 - 7-7; differentiating we have 


dr 2 
dt 


2r 


d7 

dt 
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dividing through by 2 and differentiating again we obtain 


1 d 2 r 2 

2 dt^ 



( 79 ) 


The energy integral is once again written 



= 2V - K 


where V can be written as 


v = 




(80) 


i.e., v = g hu. Upon further inspection, in the above energy integral, dp/dt remains finite and the 
last two terms in the bracket of Equation 80 are finite if r - 0 as t - 1 , . Under these conditions, 
the energy integral and Equation 80 yield 


dr V 2f ( m o +m i) 

Jt) ? L . (81) 


Lj = finite as t - . We note that dr/dt behaves roughly as the inverse root of r. Dot multiplying 

the first equation in 35 by r, we have 


_ d 2 r _ f ( m o +m i) 

r ‘ dt 2 r " L2 (82) 


L 2 - finite as t -*t r Adding Equations 81 and 82 and making use of 79 we have 


1 d 2 r 2 f ( m o +m i) 

2 dt 2 r ' L (83) 


where L' = L x +L 2 . It is clear from 83 that as r -o (double collision) d 2 r 2 /dt 2 -® and 


d 2 r 2 
dt 2 


> 0 


(84) 
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This tells us that dr 2 /dt is increasing and we have two subcases 


\ dr 2 

a) “dt" < 0 

b) > 0 (85) 


Case b) insures that r 2 is increasing so that r does not approach zero. Since we are considering 
collision we reject this case and consider only case a. In this case r 2 decreases monotonically. 

From the vector indentity a 2 b 2 = (a • b) 2 + (ax b) 2 we let a - r, b - dr/dt and write 


-2 ( d *\ 2 - (- d7 V dlV 

r \3t / - V r ' dtj + \ r * dt) 


or 


72 (ar) = r 2 (ar) + ( 7x a?) 


( 86 ) 


From 81 we see that as t - tj 


/dr\ 2 . v 

r ^dtJ ’ +2f ( m o +m i) 


and 


r 


2 



(87) 


Now from Equations 86 we see that we have the sum of two positive quantities equal to zero and 
hence they are each identically zero. Therefore 


r dt " 0 
- dr* 

r x dt ■* 0 (88) 


as t - 1 , . 
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We will now show that 


_d 

dt 


( Fx 3 f) - 0 


as t — t , . From Equation 35 


d 2 

dt 2 


f = 7, * *ft) 


We cross each side of this equation by r so that 


d 2 ? 
dt 2 


f m 2 r x p 


U - — 

1 r 3 r 3 
Vo r i 


and 


f m 2 r x p 


U - — 

W r i 3 , 


A U 

dt \ r * dt j 

Since a 3 -b 3 = (a -b) (a 2 + ab + b 2 ) we write 

it (r* 5t) = !■»,?« *77) 


or 


*M) * 


taking absolute values 


_d 

(- d7 ) 

dt 

V x dt) 


- f m 2 | r x p\ ~ 


+ — — + — 

2 r- r 2 r 3 


(89) 
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Figure 5 


In the time interval t 1 -S Q <t<t 1 , /o<r 0 +gm 1 <r 0 + r where r - min (r 0 , r x , r 2 ) also since r < r 0 

P < 2r 0 (90) 

in the neighborhood of collision and 

j r i - r o | < r ( 91 ) 

I r x p\ < tp (92) 

from 76 r 0 >b, r x > b and 


11 11 

r o b r i b 


(93) 


Substituting the inequalities 90, 91, 92, 93 into Equation 89 



6f m 2 r 2 

< 

b 3 


(94) 


We note that from 88 and 94 both rx df/dt and its derivative go to zero as r goes to zero. We now 
take two instants of time ( t " , t ' ) such that t 1 -8 0 <t"<t'<t 1 and 
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We see that 


If- dr V 

d7 V f‘ 

d 

l~ dr \ 


|l rx at j 

_ v rxd 0 < [„ 

dt 

V* &) 

dt 


from 94 



_d /_ d?\ f* 6f m 2 r2 

dt [ r * dt) dt < | - b 3 dt 


6f m 2 r 2 


We see that 



6f m 2 

< 


r 2 (t * - t") 


now as t ' - 1 , (collision) 



0 


and 



6fm 2 r 2 (tj ~ t ") 
b 3 


(95) 


Here we see that r x dr/dt not only goes to zero as r goes to zero, but it does it very rapidly since 
as r — o, t" - tj the numerator in 95 goes to zero rapidly. 

The unit vector r as 7 - o will reveal the type of motion near collision. In an attempt to in- 
vestigate this unit vector we write. 


-r x 




9 dr _ dr dr _ dr 

r dt ~ rr dt _ r at ~ r dt _ dr 
r 3 r 2 dt 


Therefore we write 


dr 1 /_ dr\ _ 

dt r 3 dtjx r (96) 

now from 95 


dr 


1 ( - d? \ 

6fm 2 (tj ~t) 


dt 

< j 

r 3 \ r * dt y r 

< 

b 3 

(97) 
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We introduce two moments of time t' and t" such that - 8 0 <t" < t ' <t x now 


£ 


dr 

dt dt 


and by 97 


f' - f" < 


£ si * < f 


* 6f m 2 (tj ” t) 


b 3 


dt 


( 98 ) 


We now evaluate the integral 


| (tj-t)dt = | t,dt-£ tdt = tj(t'-t")-(v - V) = t i( t '' t 'Vi( t '" t "x t ' +t ") ' 


from 98 


6f m 2 r i "I 

f'-f'l < “il - (t'-t") hi ' 2 (t' +t ")J 


b 3 


(99) 


as t ' , t " - t , , | r " - r ' | - 0 and by the Cauchy criterion, f approaches a limit. K we denote this 
limit by # we have 


1 im ? - < ? , l?l~l 

t-t. 


( 100 ) 


Because the unit vector is constant near collision we conclude that the motion is smooth and 
dismiss the possibility of m 0 andm, spiraling into collision with each other. In addition, (-t) can 
be substituted into Equation 35 without any effect. This suggests that the motion is symmetrical 
about the t axis. The two above conditions suggest a smooth and symmetric motion before and 
after collision, see Figure 6 



Figure 6 
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Since 


and 


From 87 


so that 


and as t - t x 


which leaves 


or 


from 101 


r r • r 

dr dr ^ dr 

dt ' dF r + r dF 



( 101 ) 


r 



2 f ( m o +m i) 


( 102 ) 


r 




2 f ( m o +m i) 


(103) 


we see from 97 that 


r 


3 



- 0 


r 



2 f ( m o +m i) 


(104) 


^ dl - ' / 2f ( m o +m i) 


(105) 


dr r dr ^ dr 

dF * T fft ' r 4 r 1* dF (106) 
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and since the unit vector r is constant near collision the last term in 106 is zero. From 100, 105 
and 106 we have 


/F af ^ ? /2f(m 0 +m,) 


(107) 


and 


S ~ l/2f(m 0 + ">,) 


(108) 


From 107, we can gain insight as to the behavior of the velocity vector near collision. 


BEHAVIOR OF "r" NEAR DOUBLE COLLISION AND INTRODUCTION OF REGULARIZATION 
VARIABLE "u" 

In order to investigate the asymptotic behavior of r in the neighborhood of collision, we write 
down 105 


J ' /2 f(m 0 +m,) 


which is equivalent to 


2 dr 3/2 

3 dt 


/2 f(m 0 + m, ) - 0 


(109) 


we see that 


2 dr 3/2 

3 dt + 


/2 f(m 0 + m,) 


< e 


( 110 ) 


if 1 1 - 1 1 1 < 7 ) where e and v are arbitrary constants and e - o as t 1 x now 


3 r3/2 " /2 f(m 0 +m,) (t , “ t ) 



2 dr 3/2 

3 dt 


+ 



dt 


Since r._. = o: also 
1 1 


2 

3 


/ 2 f(m 0 +m,) (t, 



2 dr 3 / 2 

3 dt + 


|^2f(m 0 +tn,) 


dt 
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from Equation 110 


2 

3 


r 3/2 


- f <■"„ * m, ) (t , -t) 


< <(>,-•) 


if 1 1 t - t ) < v and 


2r 3/2 

3 y' 2f K +m i)( t i ■ t) 


€ 

/ 2 f(m 0 +mi) 


if t , - t <7j as t-t,, e(t) goes to zero and 


2r 3 ^ 

3 Z 2 f (m 0 +mj) (t, -t) 


( 111 ) 


if we define 


then 


t - 1 , and 


2 1 _ _ 1 _ 
3 ]/ 2 l( m o +m i) A 3/2 


r 3/2 

A 3/2 (t, -t) 


1 


( 112 ) 


r-A( tl -t ) 2/3 ( 113 ) 


From 113 we note the asymptotic behavior of r, as a function of t, near collision. From 105 and 
113 we write 


dr jg£K + m i) 

dt ' " 
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and if we define 


B 


/2 f (m 0 + m7) 

7k 


we have 


dr 

dt 




From 115 the velocity, as a function of time, is evident in the neighborhood of collision, 
we note that as t - 1 , the velocity becomes infinite. 

We shall show that r (see 113) can be expanded into the following series 


r 


»2 ( l l “ t Y /3 + a 3 ( l 1 


t \ 3/3 + 


(*i 


- + 


if we let (tj - t ) 1/3 - u 


r - a 2 u 2 T a 3 u 3 + a 4 u 4 + * * • 


and 


du 

dt 


1 

3 


(*. _t ) 


-2/3 


or 


du j. 

dt “ r 


where 


u 


r 


dt 

r 


(114) 

(115) 
Once again 


(116) 


(117) 
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now 


du 


dt 


(118) 


from 113, 115 and 117 


dr 

du 


dr dt r 
dt dd ~ 


and 


dr 

du 


Yr 


(119) 


The regularization variable u , as defined by 117, is referred to as the pseudo-time. This new 
pseudo-time is seen to remove the singularity in the velocity, for now, the new velocity is pro- 
portional to the square root of r and as r-o, for collision, the velocity does not go to infinity. The 
velocity is now an analytic function. 

In an effort to investigate the convergence of the Taylor series solution, we state the Cauchy - 
Picard Theorem without proof. 


CAUCHY-PICARD THEOREM 

Theorem 3: Let Q, ( , q 2 • • * q n ) i = 1, 2, • * * n be analytic functions which do not contain t 
explicitly and which are developable into Taylor series in the powers of differences of q - q. and 
these series are convergent if; 


a) 


k -<*i\ < k 


Then there exists positive and finite quantities Q/ such that when a) is satisfied 


b) 


Qj (*>1 '• 'In) < Qj' 


j = 1,2, • • • n 


(120) 


Under conditions a) and b) the system of differential equations 


dq ; 

St = Qj ('ll ••• O (121) 


admits one and only one analytic solution such that q ; goes to a finite limit q. when t -t. 
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Condition a) assures the functions are analytic within the radius q/ . In this solution, the un- 
knowns ( q ; ) are developable into Taylor series in powers of t - t which are convergent for 1 1 - 1 1 -T 
where T' is minimum 


q t 52_ 5n_ 

Q,' ’ Q 2 <V 


(122) 


At a time t = t let the components of q and q be (x. , y . , z ; )^ and (x 4 , y k , z. ) respectively and 
similarly let the components of q and q be ( x t , y. , z.) and Xj, yi , z t . Condition a) assures that 


x i - x i I • |y i "Vi ■ z , ~ z i | < k o 


x i - x i | • |Vi -Vi I ■ I*i - z i | < k o' 


(123) 


where k 0 and k 0 ' are minimum radii of convergence. From 28’ 


2 = 


2 = 


(“a ~'i) + [( ? 2 _7 a) ■ ( 7 1 ' 7 i)]| 2 

|? 2 ' ? l| 2 + 2 (V ? l)[(V ? 2 HV ? l)] + [(V ? a)-(V ? l)]' 


(124) 


We let 


Eo = ( ? 2 [( 7 2 ' (^1 ~ T l)] + [(V^MV^l)] 


(125) 


so that 


2 > 


2 _ 


(126) 


from 125 


Po 






We know from 123 



< 



+ 


y 2 ~y 2 



< 3k 


o 


(127) 


(128) 
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and similarly 


r i” r i < 3k 0 


We see that 


|7 2 -T 2 | = y/(x 2 -x 2 ) 2 + (y 2 -y 2 ) 2 + (z 2 -z 2 ) 2 < fekf = V3 k„ 

likewise 


|?1 -7,| </3k 0 , |r 0 -?o| </3k 0 

Using 128, 128', 129, 129' and 127 we find 


i 12F n k f 


12 k 0 2 


and from 126 


r o 2 


- 12 ( 


r k + k ^ 
r o K 0 K 0 } 


r 0 2 must remain positive so that 


12 k 0 2 + 127 0 k 0 - 7 0 2 < 0 

Since k 0 must be greater than zero we have 

(4 /3 - 6 ) _ _ 7 o _ 7 o 

0 < k o < T2 r o " 4^3 +6 " 12.928 •" 


(128') 


(129) 


(129’) 


(130) 


(131) 


The inequality is strengthened if the denominator is set equal to 14 



and similarly 


r 


l 


> 14 k Q 


(132) 


from 124 and 125 


r 


2 

0 


r o 2 


+ p 


0 


or 


r 


o 



and 



(133) 


Here we see that l/r 0 , which appears in the equations of motion, remains analytic for r 0 > 14 k 0 . 
It must be remembered that x ; , y . , z . were chosen arbitrarily and that the radius of convergence 
k 0 depends upon these initial conditions. As r 0 is taken smaller and smaller, the radius of con 
vergence decreases. Rewriting 130 


r q 2 > r 0 2 - (l2 7 0 k 0 + 12k 0 2 ) 


using 132 


/ _ r o r o 2 

r 0 2 > r Q 2 * \12 r Q 14 + 12 195 


or 


2 _ 


12 r J 


[*] 


4 7 2 
H r o 
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so that 


and similarly 


r >^c 
r o - 7 


2 r„ 


2 r„ 


so that 


J_ < _X- 

r o 2r o 


J. < 


7 

2 7„ 


r 2 2r 0 


(134) 


These three inequalities hold for the given initial conditions 123. Using these same initial con- 
ditions the disturbing force will be investigated. We write; 


(*2"X 2 ) ' ( X 1 _X l) + ( X 2 _X l) 


Using 123 and 132 


x 2 - Xj £ x 2 ~x 2 | + Xj “ Xj + j x 2 “ X j 


8r, 


| X 2 - x i| < k o + k o + r o - 14 + 14 + r o 


(135) 


from 134, 135 and 132 


x 2 “ Xj 

< 1^0 

7 y< i 

r 0 3 

- 7 I 

r o/ 4 k 0 2 


similarly 


X 1 _x o 


X 2 -x 0 

r 3 
r 2 


r, 3 


4 k 0 2 


(136) 


Now from 129, 129' we see that condition 120 a) of the Cauchy-Picard Theorem is satisfied. 
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The disturbing potential 43 is written 


u 


/ m 0 m l mj m 2 m Q m 2 \ 

f br + ~ + ^r) 


now 


1 _au 

mj (?Xj 



and from 136 


1 iH 

irij ^Xj 


f ( m 0 +m i 


4 k 0 2 ) 


f(m 0 + m , +m 2 ) 


4k 0 2 


In a similar way, it can be shown that 


i au 


i au 



| m i 


J_ jju JM_ 
m , dz i K 4k„ 2 


(137) 


From 27, 40 and 44 


m q m j m ^ 

V -- iT~ v 


(138) 


and from 43 and the instant t - t 


u 


r m o m i 

f hr 


m 0 m 2 m l m 2 ~| 

+ J 


and using 132 


r m o 

u < f 


*i +m o 


n 2 + mi 


14k n 
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so that from 138 


n o m i m 2 _ f (m l m 2 + m 2 m 0 + nij m 0 ) 

M V < I4k7 


(139) 


Two useful inequalities will now be derived. It is evident that 


( m l ~ m 2 ) 2 + ( m 2 ' m o) 2 + ( m l ~ m o) 2 > 0 


or 


2m 0 2 + 2m 2 + 2m 2 2 - 2 (mij m 2 + m 2 m 0 + nij m Q ) > 0 


So that 


m Q 2 + rrij 2 + m 2 2 > 01^2 + m 2 m 0 + m x m Q 


adding 2(m t m 2 +m 2 m 0 + m j m Q ) to both sides of this inequality; (m 0 +m 1 + m 2 ) 2 


3 (n 


n o ) 


or 


M 2 > 


3 (ra, 


m 2 + m 2 m Q m j 


m o) 


(140) 


Now if 


4m. m j < M 2 


(141) 


then 


4m 0 mj < (m 0 + m x + m 2 ) 2 

4m 0 m t < ( m 0 + m i ) 2 + 2rn 2 ( m o + m i ) + m 2 
0 < (m 0 - m, ) 2 + 2m 2 (m Q +H1, ) + m 2 2 


43 



Since the final step is correct, then the steps can be reversed and assumption 141 is correct. 
139 and 140 


from 141 


m o m i m 2 - „ fM 2 
M V - 42k„ 


m. m. 

i j 

M 


K 



where K is the energy integral constant. The energy integral can be written 

m 0 m. m 2 2m 0 nij m 2 

— M— K + — M— V 


so that 



m 0 nij m 2 K 


2m 0 m 1 m 2 V 

M 

+ 

M 


Since this is true at any moment of time 


m q m | m 2 K 


2m 0 nij m 2 V 

M 

+ 

M 


from 142 and 143 


m Q M 


+ 


fM 2 

21k 0 


or 



fM 2 

21m 0 k 0 


From 


(142) 


(143) 


(144) 
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If we let m = minimum (m 0 , m,, m 2 ) , 144 can be written 



fM 2 

21mk 0 


or 




fM 2 

21mk 0 


It is obvious that a similar procedure will give the same results for x. , y. and z i 


so that 


x i ■ y 


^ | ■ |*i| < /t Ik| 


fM 2 

+ 21mk. 


(145) 


Since from 123 


"*i < k o' 


(146) 


We let 


k 


0 



(147) 


and 145, 146, 147 yield 



+ 



< 2k. 


(148) 


We have k 0 ' for the Cauchy-Picard condition 123. Note that k 0 ' is the radius of convergence 
for a series, in powers of differences in velocities, which represents the right side of our dif- 
ferential equations. One should realize that k 0 ' can be arbitrarily chosen. The Taylor series 
representations of the right side of our differential equations, in powers of the differences in co- 
ordinates and powers of the differences in velocities, have been investigated. The radii of con- 
vergence were also determined. Each power series representation is valid up to the singularity 
at the point of collision. A new independent variable will be defined such that the singularity at the 
point of collision is removed. 
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RADIUS OF CONVERGENCE FOR TIME SERIES WHEN MUTUAL DISTANCES ARE ALL 
LARGER THAN 14k 0 


Our final Taylor series solution for the right hand side of our equations of motion, is ex- 
pressed in powers of t - t which is convergent for t - t IT' where T ' is minimum of q,'/Qi' , 
q */ Q . . . q n '/ Q n ' . (see 121) In order to investigate the radius of convergence for the solution 

series, we take two representative equations of motion: 


dt ' *o 

di o _ .. 

dt ~ x o 

from 123, |x 0 -x 0 |<k 0 and from 148 |i 0 | <2k 0 ' so that 

q° _ ^0 

q7 = 2k o 

from 123 |i 0 -x 0 | <k 0 ' and from 137 |x 0 | < fM/4k 0 2 so that 

q/ 4k 0 ' k 0 2 
q; = fM 

We know that T' is the minimum of 150 and 151. Now 


4k 0 2 k Q ' k^ 
M ” 2k 0 ' 


k 0 

2Mkg' 


(8kg kg' 2 -M) 


(149) 


(150) 


(151) 


and from 147 


4kg 2 kg' kg 

M 2kg' 


k o r 8M 

2kg' l21m 


+ 2k„ 



(152) 


Since the largest m can be is M/3 , we substitute this into 152 and find 152 to be positive, so that 
k 0 /2kg' is the minimum of 150 and 151. We see that the series solution, in powers of t - t, is 
convergent for 


< T' 




4M 2 


r 2 ^ 14k 0 


21mk 0 


+ M|K| 


(153) 
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Motion is regular in the given interval, if the coordinates of the bodies are analytic functions of the 
time in this interval. 


This analysis has been done under the restriction that all the mutual distances are greater 
than I4k 0 . 


THE EQUATIONS OF MOTION USING NEW INDEPENDENT VARIABLE "u" 

If the scalar components of r (25) and ~p (26) are designated (x, y, z) and (£, 77 , £) respectively, 
equations 35 are written 



(154) 


The definitions of (X, Y, z, =, H, Z) should be evident from the above equations. From equations 41 
we write in scalar form. 


and from 45 we have 


g(xy -yx) + h(^rj-T^) - gh c„ 

g(yz - zy) + h(? 7 C - Jtj) = ghc, 
g(zx - xz) + h(££ - ) = gh c 2 


g(x 2 + y 2 + z 2 ) + h(^ 2 +r) 2 + £ 2 ) - 2V-K 


(155) 


(156) 
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Now from 117, 118 we remember 


dt - r du 


r 


dt 

~F 


t 


f 


r du 


With these definitions in mind we define different derivitives of a function w . Let 


dw . , dw 

dt = w and dl7 - w 


so that 


w 





w' 


and 


w ' - rw , w* = r 2 w + r rw 

Now from 160, x" = r 2 x + r rx and from 159 x“ = r 2 x + r' x'/r . Similar equations exist for 
z" so that from 154 


dx ' __ 

"SI = x 



f(m 0 +m,)x 


fr 2 X 


dy ' _ 

“du - y 



f ( m o +m i)y 
r 


+ 


fr 2 Y 


dz ' 
du 



f (m 0 +m 1 )z 
r 


+ 


fr 2 Z 


using Equations 156 and 159 we have 

g( x '2 +y ' 2 +2' 2 ) + hr 2 (^ 2 + ^ 2 + L 2 ) = r 2 (2V - K) 


(157) 


(158) 


(159) 


(160) 

and 


(161) 


(162) 
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From 160 


(163) 


using the energy integral (45) 


_ hr (dp \ 2 2 Vr K 

r r - ■ g Vdt ) + T " 1 r 


or using 27 


r2 = 


g V? 


' g 


2 ( m o +m i) f (? 


m 2 \ 
m l r J 


(164) 


Now r 2 r will be evaluated. We assert that 


r 2 r - - f(m Q +mj) + rf(xX+yY + zZ) 


(165) 


since from 154 


r 2 r - - f (m 0 + nrij ) + r 


xx + 


f ( m o +m l) , .. f ( m O +m l) , .. f ( m O +m l) 


+ yy + 


+ zz + 


f (m 0 + nij] 

= - f (m Q +11^) + (x 2 +y 2 + z 2 ) + r(xx + yy + zz) 


using 164, 165 to evaluate 163 we may write 


dr ' 


Uu - f ( m o +m i ) + rL 


(166) 


where 


L = 


2fm (m 0 +m 1 ) 2f m 2 (m 0 + m, ) h K 

f ( x X + yY + zZ) + ^ + gU 2 + V* + i 2 ) ~ ¥ 


(167) 
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We now let 


a 



f(m 0 +m,)x 
r 


n 



f( m o +m,)y 
_ _ 


7 



f (m 0 + z 
r 


(168) 


differentiating a with respect to u 


da 

dd 




i 

xr 

T~ 


Substituting x' (161) and r" (166) we have 


da 

du 



f(m 0 + mj) 

; xr'+frr 



f(m 0 +m,) _ f (m Q +m,) 


x + 


after canceling appropriate terms we find that 


a' = fXrr' + Lx' 


and similarly 


/S' 

- fY rr ' 

+ Ly' 

y' 

- fZrr' 

+ Lz' 


(169) 


From 160 and 159 


r 2 x + rf X - r 2 X + -JT x' 


(170) 


also 


a i f r 2 X — x ' - 


f(m 0 + m,)x f( m 0 +m i) x 

+ r 2 x + ; 


x' + r 2 x 
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so that 


dx' 

tkt x ' " ^ + fr 2 x 


and similarly 





dy ' 

dil ^ y" 11 P + fr 2 Y 






dz ' 

du " z “ y + fr 2 z 


(171) 

Collecting equations, we have 





dr 

HH 

r ' 

dr ' 
du 

r" - f (m 0 + m, ) + rL 

dt 

du 

r 

dx 

du 

z x ' 

dx ' 
du 

a + f r 2 X 

da 

du 

fX r r ' + Lx ' 

by 

du 

= y ' 

du 

[1 + f r 2 Y 

H 

fY rr ' + Ly # 

dz 

du 

z ' 

dz ' 
du 

y + f r 2 Z 

d> 

Hu " 

fZ rr ' + Lz ' 

df 

du 

- r f 

d 7] 

du 

r n 

d£ 

dU - 

rC 

df 

du 

= fr H 

d r] 

du " 

fr J< 

dl, 

du = 

fr S (172) 

For a given set of initial conditions 0 
from 168. In addition since 

z, x', y', z ' , f , T) , z' , T 

1 f ) a > ; 

, y can be determined 


x 2 4 y 2 + z 2 


xx ' + yy ' + zz' 


(173) 


r and r' are determined from the initial conditions. 

In order to determine L (167), we must find K from 162. It is necessary to divide by r 2 , which 
will go to zero near collision. If the initial conditions are chosen near collision, then numerical dif- 
ficulties arise. Given the initial conditions, the eighteen equations of motion (172) are now soluable. 

Our new pseudo-time variable "u” (157) must go to collision simultaneously with "t". That is, 
from its definition, ”u" appears to go to infinity. In the following exposition, which differs from 


51 



Sundman's, "u" will be shown to be limited. From 105 


; + o* o -1 


( 174 ) 


where !<! *0 as r-r 1 so that 


dt 


M*o tm i) 


( 1 + e ) dr 


(175) 


from 157 



(176) 


where the first integral is the interval near collision and the second is arbitrary since t 0 is 
arbitrary. From 175, 176 




dr 

(1 + o /F 


1 p dr 

l /2f(m 0 +m t ) J r /? 


1 f r £ dr 

V , 2 f(m 0 “+m,)J r /F ( 17 7) 


evaluating the first integral we have 


rf= — T 2 (if^-yF) 
/ 2f ( m o +m i) 

which is finite as r- r t the second integral 


1 

f € dr 

r 

Ul dr 

/2f(m 0 + m,) . 

l rr 

/2f(m 0 + m,) J r 

Tr 


there exists a number "V such that I f I < v since |e| - 0. With this in mind 

i r ui dr < i r zl^i 

y/ 2 f (m 0 + m, ) J r /F ^2f(m 0 + J r /F 


(178) 


(179) 


( 180 ) 
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The first integral in 177 is finite and the second integral, from 179 and 180 is bounded, so that "u" 
is bounded and cannot go to infinity as t -♦ . 

The question arises whether the equations of motion 29, 35 or 172 best describe the motion 
near collision. The distance r doesn T t enter into the denominator in the right hand members of 
equations 172 and this system will be generally preferable when r is small with respect to the 
other two distances. 

Corresponding to our eighteen equations 172 we have 


r , r ' , t , x, x ' , 


^ y, y', A z, z \ y, t?, <5, y, l 


(181) 


as our eighteen unknowns. If the distances r 0 , r x , r 2 are greater than zero for t - t 0 , the variables 
in 29 and consequently the unknowns 181 are developable in powers of t - 1 0 if 1 1 - 1 0 1 is suffi- 
ciently small. This follows from 117. Since r is not zero for t - t 0 , 1/r is analytic and can be 
expressed as a power series in t -t 0 . Integrating this power series from t to t 0 , u is expressed 
analytically as a power series in t - t 0 . By rewriting 117 as 


t 


f 


r du 


it is evident that the variable t is developable in powers of u when j u| is smaller than a certain 
value. Upon substituting for t , this series in the series expression for the unknowns 181, we 
obtain functions of the variable u which verify the equations 172. Furthermore, upon substituting 
the series expression for u into the solution of the system, we get the solution of the system 29 
from which we started. The equations of motion 172, are equivalent to equation 35 from which they 
were derived, and equations 35 were in turn derived from 29. The initial conditions for the system 
172, are 


x, y, z, x\ y\ z\ g, y, g\ y\ V 


(182) 


and are regulated by 162, 168 and 173. By recombination of the equations 172 the system is seen 
to have the following integrals, 


r' ( m o im l) x 


-t- 

o 

6 

1 


1 

►i 

IT 

o 

3 

N 

' ~ T x r 

:: a r 

r y 4 

r - h r 

r z + 


which are seen to be zero by 168. It remains for us to see that "r", which is amoung the unknowns 
of the system 172, as well as the constant K, which enters into the same system, both have the 
same significance as in 154 and 156. 

Rewriting Equation 173, 


rr ' xx ' + yy ' + zz ' 
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we see that this equation exists for all values of t or u and upon integrating 


r 2 - (x 2 +y 2 +z 2 ) 


is equal to a constant of integration. This constant is seen to be zero so that the unknown M r" 
satisfies Equation 32. Upon introducing "t” into 162 we get 45, from which one concludes that the 
constant of kinetic energy K is the same for both systems. 

We would now like to see the unknowns 181 developed into a series in powers of u - u x and to 
determine the lower limit of the radii of convergence of these developments. To do this, we must 
find the upper limits for the unknowns when t - as well as the absolute values of the second 
members of the Equations 172. From these values, we then can form the ratios q//Q/, q 2 '/Q 2 ' 9 * * 
needed for the Cauchy-Picard Theorem. 

Suppose now that the motion is regular in the interval 0 < t < t t where t x designates the mo- 
ment of collision. Recall the quantities K : 


r i r 2 , r 0 , r p r \ x, y,z, x y ' , z ' 


£,17. V . CL, P, y 


whose values, as t approaches are denoted 


(r)o- ( r o)o- ( r i)o- ( r 'V x 0’ Vo- z c x o ’ Vo * z o 
P O’ =0’ Po' -0’ =0 ’ ^0 ’ ^0 ’ a 0’ ^0' 


We assume furthermore that 


( r ) D 



(183) 


and 


p, - 14k, 


(184) 


where k, designates a positive constant whose value will be determined later. 

We would like to have the unknowns 181 expressed as power series inu-u, and to determine 
a lower limit of the radii of convergence of these developments. This will be the goal of the 
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following exposition. As we approach collision 


( r 0 ) 0 ’ ( r l)o > f> 1 ~ ( r >0 


(185) 


from 183 and 184 


( r o ) 0 ’ ( r 1 ) 0 > ^^1 


2 


or 


27k j 

( r o ) o 1 ( r l ) o > 2 


(186) 


also 


fM fM < 2fM ( m o + m i) 

m o ( r o)o + m i ( r i)o 27k 1 m 0 m 1 (187) 

we let 

m 2 f(m 0 +m 1 ) 2 | K j 

27m () m 1 k 1 + 4g (188) 


and from 162 


(2V -K) 


hi 

g 




using 80 we have 


X 


' 2 


■+ Z 


2r 2 fM f 1 + 1 

g L m o r o + m i r i 




2 ) 


incorporating 187 

, ,, , 2 . 2r 2 2Mf ( m o 4 m i) 2rfM r 2 K 

x ' 1 y 4 z ' g 27m 0 m, kj 4 gm 2 ' g 
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using 188 and the definition of g (27) 


x' 2 1 y' 2 4 z' 2 < 4(r) 0 2 A, + 2(r) 0 f (m 0 Hn, ) 


or 


x o ' 2 + y 0 ' 2 + z o ' 2 < 2( r >0 [ f ( m o +m i) + 2 A i ( r )o] 


it follows that 


x 0 ' | - 1 y o' I • | z o' | and | r o' 



+ nij ] + 2 Aj 



from 168 and the above inequality 


K|. \f* 0 \- and l y o| < 3 f ( m o 4 m i) + 4 A i(o 0 


using 183 


x o' | - | y o' t ■ I z o I and 




+ m,) + A, k,_ 


(189) 


and 


Ki- |/ 3 0 i and |?o| < 3 f ( m 0 +m l) + 2 A i k i 


(190) 


also from 183 we have 


x o | 4 y 0 



(191) 


If we now consider the case where ( r ) t r 0 for t - t , 


x , " Yi r z i " < r >, 


(192) 
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by using 118 with 105 


y/ z/ r j' 0 


from 107 (in scalar form) we have 


2f(m 0 1 m 


' 2f ( m o 1 IHj 


(194) 


where . , , are components of the unit vector /• and are further defined by 


using 160, 168 and 195 


a* m. + m, 


Pi - + v 


/. ■= ( m n + m . 


The values of 4, 1? ■ l are seen to satisfy 190. 

We now turn to look at the upper limits of the absolute value of the second member of Equa- 
tions 172, and state that the unknowns verify the conditions: 

: r - ( r ) J ! . x-Xjj. ! y - y ! ! . ! Z - z J j , I f - - J ! 4 | - - J | 


! r ' “ ( r ' )i : • : x ' ' x i' i ’ y ' " y/ j 


I z' - Zj' : ■ k ' 


I *> / j! U 


The constants k ' , k" , k and / ' are finite, positive quantities and we will determine their values 
later. 


Using the last relation in 28 r 

in scalar form 



r o 2 

- (^-/j.x) 2 + (v ~ f^y ) 

2 + 


which is equivalent to 




r o 2 - /'. 2 + 2 s, [ 

(?-?,) -x,) -MX,] 

M 

X 

i 

X 

1 

1 

1 

+ 2n, 

[(77-77,) -m(y -y,) -My,_ 

l + l 

[(h~ T h) -m(y -y,) ~w l ] 

+ 24 , 

[(4-4l) -m(z-z,) -MZ,] 

+ 

[U“£l) -m(z-Zj) -mz,] 


and we write 


r 


2 

0 


Pi 2 + P, 


where Pj is a polynomial in 

x - Xj , y - y ! , z - £ “ ^i’ x i» ^1’ z i 


(198) 


Since 


and 


\*ib hi I 

|C,| ip, (199) 


We see that by 191 and 197 


P,| < 12/vjk, + 12k, 2 


( 200 ) 
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and from 184 


k i 



so that 


, , 45 

P i I < AO 


49 P l 


(201) 


We now see from 198, 184 and 201 that 1/ r 0 is developable into a series, of powers in x - x 
y " y i > z " z i> > V ~V 1 and £ - . From 198 


i 2 _ 


and using 201 


2 

> 7 P i 


( 202 ) 


Similarly, if we start with 


r j 2 = (j+^x) 2 + (rj + \y ) 2 + (5+\z) 2 


we can show 


2 

> 7 


(203) 


Using 184, 202 and 203 become 


> 4k, 


(204) 


Here it is evident that both l/r 0 and 1/r, are analytic in the interval before collision. We state 
that 




V 






V 


r, 

15 

r o 


r 3 
r o 






r l 3 


^ 7 3 

64k j 2 
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1 

64k, 2 


(205) 


X 


y 


z 


X 


y 


z 

r 3 
r o 


r 3 

1 o 


r 3 
r 0 


rj 3 


r i 3 


^T 3 


The first inequality is shown as follows: from 197 and 184 


whereby from 199 

15 pj 

1^1 < 14 + P\ = “TT 

and by using the above inequality and 202 

J_ < 735 
r o 3 16 Pi 


or finally by 184 

15 

64k x 2 



The other inequalities of the first line in 205 are handled analogously. The inequalities involved 
in the second line of 205, follow from 204, and 191, 197 which together show that 


i X j , I y | , 1 z | < kj < 206 ) 

Using 154 to define ! fX| , | fY| , | fz| , | f=l , 1 fK| , | fS|, remembering that the sum of absolute values 
is greater than the absolute value of the sums, and employing 205 we find, 


m 2 

Ifxl, I f y) , Ifzl < — 2 

Usl , ! fK| , I fS| < ^ 

Now using 167 to define |L| , and employing 204, 206, 207 we note that 

|l| < x, 


(207) 


(208) 
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where 


where the factor 6v 


Now from 183, 


from 208, 211 


from 207, 211 


from 207, 208, 211 


fX r r ' 4 Lx 


from 211 


m 2 r 

2k, [3* 


^ m n m 7 ' i ( V l 2 * 6v i k i' +3k i' 2 ) * IT 


(209) 


k,' + 3k,' 2 has been arbitrarily added and the velocity v, is given by 


yl i + V * i , 2 


( 210 ) 


189, 190, 197, 210 


r < k. 


U' ! . Iy'| . I z' I . I r ' | < k' + ^k, (m Q + m, + A, k, ) 


!?l, Id < v, + k" 


|aj , 1/4, \y\ < k + 3f (m 0 + m,) + 2 A, k, 


( 211 ) 


m 0 +m 1 +rL < m Q + m j + kj 


4 r 2 fX| , \0 4 r 2 fY | , I / + r 2 fZ| < k + 3f(m Q + ) + 2 \ kj + 


~T 


fY r r ' 4 Ly ' | , ] f Z r r ' + Lz ' | < ^k { ^ [k ' 4 j/kj (m Q 4 + \ 1 k^j 


■■f , | r // ! , | r f, | < k i ( v i + k " ) 
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from 207 and 211 


■fHl, IrfKI. |rfZ| < 


Returning to the system 172, we conclude that the quantities which correspond to the ratios 
q '/Q - ,q||o 2 ' + • • • , q„'/o n ' in the Cauchy-Picard Theorem, in the present case are; 


>jk' + ^kj (m 0 + nij + Ajkj)J 


m o +m l + K l k i 


k + + 3(m 0 + rrij ) + 2 A, kj 


r,) [ k ' + Al ( m o +m l + A i k lj] 


2 ( v i ■* k" 


Since ' can be arbitrarily fixed, we give it a value such that r'/k l is larger than the other ratios 
and we choose for k ' , k" and k the values 


^kj (m Q + m j + Aj kj ] 


1 /M 

2 Vk, 


k 1 m o + m i + \ k i 
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These choices of k' , k", k bring the above ratios to the following form 


/ k i ( m 0 1 "V \ k l) 

4(m 0 + m j + Aj kj j 

✓M- o + m i + A 1 k i) 

m 0 + k l 

Ai (m 0 + m, + A, k,) 

__ 

T + 4 ( m o + m i) + 3A, kj 


Ai ( m o + m, + A, k,) 
2Aj k x + m 2 


1 

JW 

2v > + At 



We shall designate by Q 2 ' the smallest of these quantities. Then by the Cauchy-Picard Theorem 
we know 

1) that, in the solution of the equations 172, which were deduced from 29, the quantities (K), 
are developable in series with powers of u -u 1 

2) that these series converge for 


l u " u i| * (213) 

3) that the inequalities 197 will be such that "u" verifies the inequality 213. 

We now consider the case where one of the distances (r) goes to zero as t — t x . From 50 


such that 


1 im R - R. 

t . 1 


R, = /h/', 
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and from 184 as t - t x 


so that 


As t - t j , r ♦ 0 and 


( r o ) i ( r 1 ) 1 


Using the Equations 172, and the method of undetermined coefficients the following series 
solutions are obtained. These results can also be obtained using the asymptotic behavior of dr/dt 
and dr/dt . Only the first terms are given. 


f( m 0 +m 1 ) • 


r ! + 6 f| ( u “ u i) 3 + " 


Mf (m n + m.) 

e — A (u- u \3 + 

fin 3 V U U ll 


M + 6 ^1 ( U_U l) + 


Mf (m n +01.) 

*>i 771 ' ( u_u i) ; 


f ( m o im i) • . . , 

6 M u ' u i) 3 + 


Mf ( m o +m i) . 

771 "h ( u_u i) + 


f(m 0 + m, ) 


"2 ^( u_u i) * 


f (m 0 ^mj) <#>(u-u,) + 


f K +m l) , 

y " 2 y\ u_u i) 


y' " f ( m o +m i) x( u-u i) + 


f ( m 0 +m l) , 

z • 2 7 (u -u, 


f K 4 m i) u i ) + 


2 ( u_u i) + 


f( m o + m,) (li-u,) + 
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t 


t 


(ll -II ,) 3 


I • • 


(217) 


f(m 0 i m,) 


Looking at the series 217, we can see that u -u 1 can be developed into a series of (t - tj 173 
and substituting this series in place of u - Uj in the formulas 216, we find that the quantities f , r h 
:, • • * are also developable in powers of (t - tj 173 . The quantities u, «f, r h • • • considered as 
functions of t , have singular points at t - t 1 . 

The same series can be used to describe the movement of the bodies after collision. 

From 216 we note that the ratios x/r, y/r, z/r tend toward the same limits (i.e.,0, when 
t goes to tj in an increasing or decreasing fashion. We must conclude from this, that the motion 
approaching and exiting from collision forms a cusp. Of course, investigating the motion of a 
colliding body, after collision, is only of mathematical interest and the fact that any motion at all 
can be investigated after collision points up a divorce of mathematics and physical motion. 


DETERMINING A LOWER LIMIT FOR Tf R M -DOUBLE COLLISION 

It will now be our purpose to reinvestigate R, (Equation 50) with the idea of finding its lower 
limit for the case when the constant of the area integral is not equal to zero (two body collision). 
From Equation 50, we see that, in the case of two body collision, R goes toward a finite and non 
zero limit. If we define an interval t ' < t < t" , then from 75b. 


H" < H' 
R' £ R' 


and from Equation 71 


H" 



+ KR 



2 


so that 


H' ^ R" 



+ KR 


+ 


g 


2 h 2 c 

R" 


2 


(218) 


and 


K 2 h 2 c 2 

R" 


H' ~ KR" ' H' 


I K | R' 
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and finally if R doesn't have a minimum for t = t * 


E 2 h 2 c 2 

R " - H' + I K | R' 


(219) 


The inequality 219 remains valid until t* goes through t, where dB/dt changes sign, and for 
the minimum value of R' , where dR'/dt - 0 , the equality 


H' 


R' 



+ KR' 


+ 


g 2 h 2 c 2 
R 1 


gives us 


g 2 h 2 c 
H' + |K| R 


for K > 0 or K < 0. We see that the relation 219 holds for the minimum value of R so that we can 
generalize and say, Theorem 4 : 


g 2 h 2 c 2 

R " H' + |K| r' 


( 220 ) 


is valid up to the point where R - R nax . 

IfK<0, Equation 51 shows that d 2 R 2 /dt 2 is never negative from which one ; concludes that « 
doesn’t have a maximum. In this case, R goes to infinity when t goes to infinity, and 220 therefore 

gives a lower limit to R which remains valid for all values of time. If K >°.^ R does ^ ^ticular 
maximum for finite values of t. Equation 220 will again give the lower limit forR . This particular 
uS?”r be expressed in the ‘following manner. We multiply the numerator and denominator of 
220 by R and employ Equation 71, to get 


R > 


E 2 h 2 c 2 R 


(£) 


+ 2KR 2 + g 2 h 2 c 2 


If we let R° and dR°/dt be the values of R and dR/dt for t 0 and let 


> 2 h 2 C 2 = f , 2 


( 221 ) 
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we have Theorem 5: 


R > 


f 2 R° 





( 222 ) 


if K < o and 


> no maxima 


R > 


f 2 R° 

/ HR °\ 2 

\R° dpj + 2KR 02 + f* 


J 


(223) 


if K > 0 and R has no maximum for finite values of t . 

It remains for us to find a lower limit of R when K > 0 and R has a maximum for a finite value 
of time (i.e., M S" expression 289). We will find this lower limit on R, for all values of time, such 
that it only depends upon andK. We assume that R has a maximum R ' for t = t ' so that we 
consider the case where 


K > o 


(224) 


and 


dR ' 
dt 


0 


(225) 


Now from 71 and 225 at t = t ' 


H' 


KR' 



(226) 


Since R' is a maximum, there exists in the neighborhood of t' , and instant t" such that the deriva- 
tive d^/dt does not change sign and that R <R' in the interval from t ' to t" . From 218 and 226 


R" 



+ KR" 


+ 



< KR' 



(227) 


so that 


K(R' -R") > f, 2 (ip 


f, 2 (R' -R") 
R" R' 


67 



and 


KR' 2 > f, 2 

and finally by taking the square root of both sides and multiplying by f , we have 

f t 2 

TT < 


(228) 


(229) 


now from 226, 229 we have 


H' < f, VR + KR' 

Starting with 223, and using 227, 226 and 230 respectively we see that 


R > 




f, 2 


„ /dR"\ 2 

* (-ar J 


f. 2 


+ 2KR" + 2KR' + 


R' 


f. 2 

H' + KR' 


f , 2 

f , /K + 2KR ' 


(230) 


or 


f 2 

R > fj /K+ 2KR' 


(231) 


Now this lower limit on R will become smaller as t goes to ±® if R' gets larger as t goes to *». 

We wish to find a positive, fixed limit which remains valid however large the maximum R' . 
With this in mind we note 156 and write 


2V - K > 0 

and from 80 


m 2 r 2 2M 


(232) 
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where f is set = 1. If we designate r ra as the smallest distance r. , we have from 232 


J_ (JL J_ J_\ > _K_ 

r n, \ m o + m i + m 2 / “ 2M (233) 


Now we note that r < q where 

m * 


2M f_l_ + _L + J_\ 

K \m 0 + nij + m 2 ) ( 234 ) 


The movement of the three bodies will be such that one of the distances r t will be the smallest. 
Since the distances r. are continuous functions of time, it is evident that each time that a certain 
distance ceases to be the smallest, it will become equal to another distance. At this instant, both 
distances are considered the minimum distance and 234 insures they are both <q. The third dis- 
tance of the isosoles triangle would be < 2 q. To simplify our following formulas we will say that 
all the distances are <qis. From 52 and the fact that r 0 , r, , and r 2 < is q we write 


R < R„ 


(235) 


where R 0 is denoted by the positive radical of the equation 


R o 2 5q2 (m 0 + m, + m 2 ) 


(236) 


We conclude that in an interval of time where R >R 0 a single distance remains constantly <q. 
We now, consider the movement in an interval of time where the inequality 


R > R„ 


(237) 


is constantly verified. We also let r 2 » r be the distance which remains smaller than q, or 


r < q 


(238) 


We will now deduce from 237 and 238 some other needed inequalities. We shall also determine 
a certain value R 0 of the quantity R which will play an important role. With this in mind, it can be 
shown that 


R „ 2 ~ (g + o- 2 h)q 2 


(239) 
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where g and h are defined in 27, 40 and 


r 2 = 4 + 


4m 0 m 1 m 0 2 + m Q m t + m 2 

, 2 ( m o + m lj + (m 0 + m,) 2 


These relations are proved by substituting 240 into 239 and winding up with 236. We can 
240 that 


a > 2 


From 50, 237, 239 and 238 respectively 


gr 


2 + hp 2 > R 0 2 = (g + o^h) q 2 > gr 2 + cr 2 hq 2 


or 


p > aq 


and from 238 




We see from 185, that r 0 >p - r, r, >p - r and from 243 and 242 


r o > ^o- — P > (o--l)q 
r, > P > (cr~ 1) q 


Rewriting Equation 50 as 


hp 2 = R 2 - gr 2 


so that 


P - 


R 

ih 


(240) 
see from 

(241) 


(242) 


(243) 


(244) 


(245) 


(246) 
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Using 239, 237, 238 


gr 2 < gq 2 


g R o . gR 2 

g + a 2 h g + cr 2 h 


(247) 


we now can see from 245, 247 that 


h p 2 > R 2 


gR 2 

g + a 2 h 



g 

g + cr 2 h 


] 


R 2 


a 2 h 
g + r 2 h 


] 


or 


P 



(248) 


where the constant 


■/ 


g + cr 2 h 


(249) 


is smaller than one. 
We now define 


R„ e 


(250) 


where we see that R 0 >R 0 . Designating by p 0 , p 0 the values of p which correspond to the values 
R 0 and R 0 of R, we see that after 246 and 248, 


P o 




(251) 


and consequently p 0 > p 0 . One must conclude from 251 that in the interval of time in which R de- 
creases from R 0 to R 0 , there exists an instant t where the inequality 

3t < 0 (252) 
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is verified. That is to say, that p is decreasing at an instant t . Letting R, p, • • • designate the 
values of R , p, • • • for an instant t , and using 252, 250 and 251 we have 


dp 

Ht < 


(253) 


R„ £ R £ R„ 


(254) 


Ro _ . 

< O — 


Vh 




(255) 


The differential equations of the movement (35) remain invariable when "t” is replaced by t . 
From this we deduce that the lower limits are independent of t and are valid before and after the 
maximum of R . We shall study the values of R after the moment t ' where R passes through a 
maximum R' . It will be necessary to divide the maxima into three classes, according to the sizes 
of the maxima R' and the minima R" which follow it. 


The first class will refer to the maxima which verify the condition 


(256) 


and after 231 we have 


R > 




f , Vk + 2KR 


0 


The second class will include the maxima for which 

R' > R 0 


and 


R" £ R 0 

and since R is in the interval of maxima R' and minima R" , 

R > R 0 


(257) 


The first two classes are seen to have lower limits. The third class satisfies the conditions that 

R' > R 0 
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and 


R" < R 0 


(258) 


The next section will be concerned with finding a lower limit for R in this third case. 

Consider a maximum of the third case. R will diminish constantly from R' to a value smaller 
than R 0 since the minimum is less than R 0 . This third case corresponds to the relations 253, 254, 
255. To find an inferior limit forR, we will seek a superior limit for H (Equation 71) when t = t . 
In order to accomplish this, we must find a limit for the absolute value of the derivitive dR/dt . 
Since from Equation 54, 


- dR _ dr _ dp 

R dt = dt + h P dt (259) 


we see that the superior limits for the expressions 


_ dr 
r St- 


and 


_ d P 

p at 


are necessary. In an effort to determine the upper limit of the former relation, we note that 



(260) 


and 



(261) 


and from 156 it follows that 


g 



< 


2V - K 


(262) 
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If R>R 0 (237) we see from 244, 241 and 234 


2M 2M 2M (J_ _1\ _ m o + m i K 

m o r o + m i r i < “ !)<1 \ m o + m i j J_ + J. + _L a ” 1 (263) 

in q ni j m ^ 


Using the definition of v (Equation 80) and setting the constant f = l for simplicity, we see 
from 263 that 


2V - K < 


2M 


(264) 


Now incorporating 264 into 262 we have 




< 


2M 


r 


or 


dr\ 2 2M 
dt J < gm 2 r 


so that 


_ dr 
r dt 



and from 234 


_ dF /2Mq~ 

r dF < y (265) 


It remains for us to find an upper limit for \p dp/dt| . To this effect we write the definition 
p 2 = f 2 + t)' 2 + i 2 and differentiate twice with respect to t and find that, 


d 2 p 
P ~dt 2 ' 


(dp \ 2 

d 2 £ d 2 rj 

d n 4 

(d £\ 2 

, (d V \ 2 


(at) 

? dt 2 v dt 2 

dt 2 

(at) 

+ (at ) 

+ lar J 
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and after 261 


p 


d 2 p 
~dt 2 


1 € & + 


d 2 rj 
dt 2 


+ 


4 


d 2 j 
dt 2 


so that by means of 154 


d 2 p 

p l^ 


> 




+ \u.M 




(<fx +rjy + £z) 


or 


p 




p 2 -/x(xg +y7? + zj) 


p p 2 + K(xg + yp + zQ 


(266) 


Now writing 28 T in scalar form 


r 0 2 = (<f“P*) 2 + (p-py) 2 + (£‘pz) 2 


we have 


r Q 2 = p 2 + p 2 r 2 - 2/x(x£ + yp + z£) 


and it follows that 


p 2 r Q 2 - p 4 + M 2 r 2 p 2 - 2p p 2 (x£ + yp + z£) 


or 


p 2 r 0 2 = [p 2 -p(x^ +y T} + zl)] 2 + M 2 [r 2 p 2 - (xf + yp + z£) 2 ] 


(267) 


now since 


|x<f + yp + zi\ < rp 
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it follows that 


p. 2 [r 2 p 2 - (x£ + yp + z£) 2 ] > 0 


and from 267 


pr„ > p 2 -p(xf +yp + z£) 


(268) 


In an analogous manner 


> p 2 +X(x£ +yp + z£) 

P r l ~ 


(269) 


Now from 266, 268, 269 


d 2 P 


> 





and 


d 2 P > 
dt 2 



Now from 244, 


d 2 P M(\ +p)(7 2 > 
dt 2 (cr-l) 2 p 2 


or 


d Ip 

dt 2 



(270) 
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where 


C 


Mcr 2 

If- Tj* 


(271) 


Before going further, we must consider separately, 

1) The case where dp/dt < 0 when t decreases from t ' to t . 

2) The case where there exists jrn instant t " between t' and t such that dp/dt = o for t = t 
and dp/dt < o between t ' and t . These are two subcases of 258. 

In the first case, between t ' and t , we see that dp/dt < 0 and from 270, 


dp fd 2 p c " 

VJ 10 


(272) 


Integrating between the limits t' and t 


J> 


d 2 p 
"dt 2 


dt + 2 


C fi c dp 

L^ dtl ° 


we have 


and finally 




2c 

p 7 



(273) 


In the second case we integrate 272 between t'" and t, and noting that dp/dt = o for t - t'" we 

have 


3t 



2c_ 

p'" 


We can now see that the inequality 273 is valid in both the first and second cases. We shall 
now return to the determination of an upper limit for | p dp/dt | . R was, by hypothesis, a maximum 
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for t = t ' so that d 2 R 2 /dt 2 <o for t = t ' and after Equation 51, we see that 


(274) 


where V' is the value of V when t = t ' . From 274, we note 2 V' - K< 2 K-K-K and from 262 


g 




(275) 


Since dR'/dt - owe find by differentiating Equation 50 

. dL.' , . . d£_' _ 

gr dt + h P dt ~ 0 


or by rearranging terms 


/dr '\ 2 hV 2 /dp '\ 2 

B Vdt ) = V dt ) 


substituting this into 275 




< K 


or 




< K 


and 


/d P '\ 2 < Kgr ' 2 


after 238, 258, 255 


dp '\ 2 < Kgq 2 Kgq 2 

•"at ' ~ hR 0 2 h 2 p 2 
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! 


from 273 


(dp) 2 < 2c Kgq 2 

' dt ' P h 2 p 2 


multiplying both sides of this inequality by p 2 and taking the square root, we have 


- dp 

P dt 




< ,/5l3- + 2c p 


and after 255 


- dp 
P dt 




2c R 0 

Yh 


Now from 259 with the aid of 265 and 276 


- dR 
* dF 


+ "j Kgq 2 + 2c R„ hVh 


where it follows that 

R 2 (ar) < + 2Kgq 2 + 4c R 0 h/h 

since R >R 0 (254) 

< T 0 + 2Keq 2 + 4c R„ hih) 
Upon recalling Equation 71, and since R > R 0 we have, 

H + KR < S 2 


(276) 


(277) 
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where 


S, = -^(^f *f, 1 *K S q'*4c6„h/h). 2K5 0 

f? = g 2 h 2 c 2 


( 278 ) 


from 75b we note H + KR-H' +KR' and from 220, 2?7 


R > 


f , 2 


This is the lower limit for the third case, 258. Since the above inequality is true for all R, one 
sees that 


R 


0 


> 


f , 2 

s 2 


Now if we note the limits for the first and second cases (256, 257) and from the definition of s 2 (278) 
we see that the inequality 


R 


f , 2 


S 2 +f, /K 


(279) 


is valid for all three cases, (256, 257, 258). The next few pages will be devoted to further rear- 
rangement of this lower limit onR. 

We shall designate by "m" the smallest of the masses m 0 , m l9 m 2 . We note right away that 


M > m 2 — m 


m 



m Q + nij > 2m 


(280) 


from 27, 40 we find 


g 


J_ < _3 

+ m 2 - 2m 


80 


h 


(281) 


— + JL <1 

m Q rrij m 


from 280 l/m > 3/M and l/m 0 + l/m 1 + i/m 2 > 9/M so that 


_9 < J. _1 JL < 3 

M “ % + m i + m 2 “ m (282) 


from 234 


q 


2MM J_ _1\ 
K \m 0 + nij + m 2 j 


and since 9/M < l/m 0 + l/n^ + l/m 2 we have 18/K<q and since l/m 0 + i/ mi + l/m 2 < 3/m we have 


18 . 6M 

K — q — Km (283) 


from 236 


~q 




+ 


l 


) 


and using 282 


3 



< 



(284) 


from 284, 



and 


_ 1 _ 

R o 
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using 283, we find 


R„ 



Now from 284, R 0 < q Vl5/m and 283, 



(285) 


(286) 


from 236, 239 


5 





g + cr 2 h 


Using 40, 240 


4 ( m o +m i) 4 m o 2 +m o m i +m i 2 

m 0 nij m 2 (m 0 + m ( j m 0 mj 


4(m 0 m 2 +mj m 2 + m 0 mj) 


m Q 2 + m 0 m i W 
+ ( m o +m i) m o m i 


so that from 249, 250, we have 


Ro 

R o 



sf-l + -1 + — ) 

3 \n> 0 m 2 / 


44 +±) 

\m 0 m, m 2 j 


m Q 2 +m 0 nij + m 2 
m o m l ( m 0 + m l) 



from 241, 271 we see that 


C < 4M 

Now from 278 and the inequalities 280, 281, 283, 285, 288. 


(287) 


(288) 


/ 8M f, 2R \ /M 

b + ^~)y 5 


2 (8 yio 


15/5) ~ 

m/m 


, and consequently 


f, 2 

S 2 + f/K > 


/ f 2 Km 2 \ /TjjT 

(8M+ 54 - Jyj + 2(8/i0 + 15 /3)M/S + fj /K m 2 
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by replacing m by M/3 in the denominator and simplifying the numerical coefficients, we have 


fi 2 


f, 2 


s 2 +f, /K 


i, 1 

64 + 9 VKM + "Jq 24 f i 


KM 


m/m 


and finally 


S 




l 

m/m 


(289) 


While considering the maxima of the third class, we have supposed that r 2 remained smaller 
than q for R>R 0 . However, since S is symmetrical with respect to the three masses m 0 , m„ m 2 , it 
is evident that the result we have obtained remains true in the cases where r 0 < q and r , < q for 
R > R 0 . In short we will have R > S in the case where K > 0, f > 0 and R has at least one maximum for 
a finite value of t . Referring back to 222 and 223 we are now ready to state an important theorem 
of Sundman. Theorem 6 : If the area integral is not zero for the three bodies one will always have 


R > J 


(290) 


where J designates the quantity 



+ f , 2 


if K 1 o and the smaller of the quantities S , and 


f 


2 

1 


R° 



+ 2K R° 2 


+ f 


2 

1 


if K> 0 . 

The case of K > 0 has two limits corresponding to the cases where R passes through a maximum 
(s) and the case where R doesn't pass through a maximum. 

The next section dealing with Sundman's exposition will treat the problem of finding a lower 
limit for the radius of convergence of the development into powers of u -Uj . 
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DETERMINATION OF A LOWER LIMIT FOR THE RADIUS OF CONVERGENCE FOR SERIES 
SOLUTIONS IN POWERS OF u - u t 

Theorem 7 : If f , > 0 , the two larger of the distances r 0 , r, , r 2 remain constantly superior to 
the quantity 


l 



(291) 


To prove the theorem we note from 52, that 


R 2 


< 



+ r 


2 

2 


If the theorem were not true, at least two of the distances r 0 , r t , r 2 would assume, at a certain 
instant, values smaller than or equal to*, (i.e., say r,, r 2 <1 ). The third distance would be less 
than the sum of the other two, and hence it would be less than or equal to 2 i. (i.e., r 0 < r, + r 2 < 2l) 
Now 


f o' W* r 2 < 6^ 


h 


2 


so that 


R 2 < | J 2 


This is a contradiction to the preceding theorem 290, and our present theorem is proved. We shall 
now fix, in a convenient manner, the constant kj . Since r 0 and r 2 are the two larger distances, 
they will be considered greater than -t (291). Noting that p > r 0 - r and p > t i - r , we conclude from 
183 and 184 that 


p > 14k t 



(292) 


Now it is clear that P > r 0 - r or p >1 -k,/2 , and if we let 4 = k,/2 + 14k, , we arrive at the first 
inequality of 292. Therefore from 291, and fixing the value of k 0 from 132, 183 we have: 


k i 



and 


, k i . lift 

k o " 28 ' 1218 


(293) 
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Under the perturbative effect of the third body, the remaining two bodies can asymptotically 
approach each other as t goes toward infinity. The question arises, as to the nature of the velocity 
( ) as the two bodies approach each other. It will be shown in the following subsection that the 
velocity v = dp/dt remains below a finite limit when r <k 1 /2 , r 0 , r, >1 . 

The proof begins by defining 


A 


( m o +m i) 

2t 


4M 


(294) 


with the assumptions being, 


r 


< 


2 


while 


> l 


(295) 


from 45, 80, 295 



I K | 


/dir\ 2 /dpV < 4M m 2 ( m o +m i) I K I m Q 2M 

g \dt / + h \dt / “ m 0 m l m 2 | 2% + 4M ^ + m 2 r 


from 27, 40, 294 



< 4gh A 


2M 


(296) 


now 


hv 2 1 4gh A + 


2M 
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and 


v 2 — 4g A + 


2M 


n o + m i 


= 2g 


m 2 (m 0 +H1,) 


2 A + 


m„ + m , 


and 


/m« mj \ 
v 2 < 2gy + 2 A ) 


multiplying both sides of 297 by r 2 and noting r < k | /2 we see that 


rv < ]/gk, (n^m, +Akj) 


Starting with 296 we can easily show that 


dx 

r dt ' 


dy 
r dt 


and 


dz 
r dt 


< ^2hr (iHqIHj + 2Ar) 


We now assume that when r < kj/2 


v 2 > D 


where 


D = 



It will be our objective to show that the assumption 299 leads to a contradiction, 
assume the inequalities 


k i 

r < T 


rv < |/gk, (m 0 m, + Ak,) 


dx dy 

r dt r dt 


(297) 


(298) 


(299) 


We therefore 
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r dF < /hk, (m 0 iHj + Ak,) 


are verified when v 2 > D. From relations 41; 300, we note 


similarly 


(fV - r^') - gc 0 ~ h 


f dy dxl 

L x Ht ~ y hi J 


- g | c o | 


g f dy dxl , | /^i , v 

h L r dt + r dt J < g | c o | +2 y~h ( m o m i + Ak i) 


l 1 ?S' + ^'l < g)c 1 |+2l/-j7(m 0 m 1 +Ak 1 ) -B 


S | c 2 + 2 Vh m o m l + Ak i 


so that 


It can be shown that 


fV - v€' I < A 


l'] < B 


v < c 


~ dpy /_ d/?y 

p • dt j + v px dt y 
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or 


p 2 (gf ) = (p ■ p ) 2 + (tv-vh 2 + ( vi-tv ) 2 + (d-d ) 2 

now since p - v and p>l- k,/2 

p 2 (a?)’ > “ A 2 - B 2 - C 2 


so that 


where 


dp 

p dT > » > 0 


'v 2 (-e- 3) -a 2 -B 2 -c 2 


is verified if v satisfies the conditions v 2 > D 


v 2 t) - A 2 - b 2 - c 2 > 0 


Employing Equations 154 we see that 


( dp\ 

d 2 p 

/dp \ 2 

, ~k p 2 -p(x£ + y-T) + *0 M p 2 + + yr) + zol 

\PSi) 

P dt 2 

vat) 

“ "L> '• r,= '1 J 


dv dp d 2 p _/d 2 1 d^J 7 dH\ dp 

dt _ 2 dt dt 2 2 \dt 2 dt 2 dt 2 / dt 


- + 71) p af + " 71) ( x ^ + w + 
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from 268, 269 


P 2 + yi7 + z£)| - P r 0 

p 2 + ^(x<f + y rj + z£)| £ p fl 


now 308 gives us 


A I dM 

dt \P dt) 


> v 2 



JL £.\ 

r i r J 


Since r 0 >p - r > p-k x /2 and v x > p- r > p- kj/2 and remembering p > I4kj (292) we have p/r 0 and 


_p_ 

r, 





28 

27 


so that by using 295 


_d 

dt 



(310) 


where 


E 


28 M 
27 Z 


(311) 


from 309 we see that dv 2 /dt will have the same sign as -p dp/dt if 





(xf + yv + zO 


(312) 


now since 


then 312 is true if 


/>3t > (x£ *yv + zb 


from 300 


xf +yT) + zU 


< rv < gk L (m 0 in, + Ak,) 


now using 305 and if we assume 


W > j/ gk, (m 0 m, + Ak,) 


(313) 


then 


p gf > W > |/gk, (m 0 m, + Ak,) > rv - (x<f + yv + zj) 


In short, we can say that if 313 is true then 312 is verified. We can go further and say that 
dv 2 /dt has the sign of the quantity -p dp/ Ax. if 


V 


2 



- A 2 - B 2 - C 2 


W 2 > gk, (m 0 m, + Ak,) 


and 


v 2 i D 


(314) 


We now let 


. it) , j/A 2 +B 2 +C 2 + gk, (m 0 m 1 + Ak,)j 


(315) 


If we choose 


G 2 = iTE, 
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and 


V 


> G„ 


(316) 


it is immediately obvious from 310 that 


d / dp\ 

di V° HtJ > E (310') 


If G 2 = r^D or 


2l - k, /A 2 +B 2 + C 2 + gk, + Ak,) 


then since it must be the maximum of the three, it follows that if G 2 = /d then D> 2E and if 


G 2 " 2-t ~ k, V A 2 + B 2 + C 2 + gk, (m 0 m, + Ak,) 

then 

l&hr] [a 2 +B 2 +C 2 + gk, (m 0 m, + Ak,)] >2E 


so that 310 T is valid for either of the 3 possible values of G 2 . 

We now try to get a contradiction to 316 to demonstrate that v always remains smaller than G 2 . 
Now there will be an instant of time t ' when v has a finite value v'(>G 2 ) and from 305, we can 
conclude that p dp/dt admits for t = t ' a finite value p' dp'/dt which verifies one or the other 
of the inequalities 


dp' 

p "gr < ~ w < 0 

dp' 

P -flt > w > 0 


Suppose that p' dp'/dt < -W'. We increase the time t after t'. Under our assumption 316, v 2 > 0 
since G 2 is, and in this first case we have assumed p dp/dt < 0. From the discussion following 313, 
we know that dv 2 /dt and -p dp/dt have the same sign so that dv 2 /dt > 0. In this first case we con- 
clude that v 2 increases positively if p dp/dt <0 (if v' > G 2 ). Since v >G 2 is assumed true we re- 
member that 305 is also verified and the desired contnidiction is reached. Now since the equations 
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of motion are invariant with respect to a change in the sign of t , we can show that the second 
inequality also leads to a contradiction for t decreasing after t' . We now conclude that the assump- 
tion 316 was incorrect and 


i- - r 2 < ~2 


We now direct our attention to look for a upper limit on G 2 which doesn't change by a per- 
mutation of m 0 , m 1 , m 2 . With this in mind, some inequalities will be developed. From 281 and 
some previous inequalities it can be shown that 


g 1 L 


m o m i < 4 


2 _ m o m i M 

h ~ m 0 + mj 4 


m o m i M 2 

m 2 m Q + rrij 4m 


M _M_ , , 
B A < 2l + 16 


m o m i m 2 , M 2 
M 1 27 


g(m 0 m 1 + Ak i) < mU 9 + 4m + 16 + ' K 


1 M K i 


4M / 1 M K i 


D < k 1 \29 + 4m + 16 |K 


A 2 + B 2 + C 2 < 2g 2 c Q 2 + Cj 2 + c 2 2 + 


V K m l +Ak l)] ( C 0 2+C l 2+C 2 2 ) 


3M 2 k j / j m 


+ m V29 + 4m + 16 |K 


Now, from the above inequalities it is found that the expressions V2E , Vd and 


2 k' A 2 + B 2 + C 2 + gk, (m 0 m, + AkJ 
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are all smaller than the quantity 



(318) 


considering the definition 315, it follows that 


G 2 < G 


(319) 


With this in mind, one notes that in order to calculate an inferior limit for the quantity Q 2 ' , 
(i.e., the smallest of the quotients following 212) we take v, = G. Thus, we now have an appropriate 
value for v, in 210 and in the quotients. It can be shown that the denominators 


4(m 0 +m, + A x kj 

m 0 + mj + kj 

m 2 

~2 + 4 ( m o +m i) + 3 A i k i 


2 \ j kj + m 2 


of the quotients are all smaller than the quantity 


4M 


5M 2 M 9 M / M , , 

■5T + m G 2 k, + 3-G^ik; +2 iKl k i 


now since 



we see that 
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is not the smallest quotient and is eliminated as a possible lower limit for the radius of convergence. 
Since 


m 


0 


+ A, k. 



or 


Ai (m 0 +m l +A l k l) 



it follows that all the quotients following 212 are greater than 


Q 


6 + 


15M 

8m 




3 

4 


|k| k, 


or 


1 


2G + 



(320) 


whichever is the smallest. The first expression is always the smaller of the two, and consequently 


Q 2 ' > Q 


(321) 


We now conclude that the development of the unknowns of the Equations 172 into powers of 
u - U j converge if u verifies the condition 


|u-Uj| < Q ( 322 ) 

Hence, we have succeeded in establishing a lower limit for the radius of convergence of the series 
solutions for Equations 172. 


INTRODUCTION OF A NEW INDEPENDENT VARIABLE "w" 

In the preceding work, we have employed in the place of t , a pseudo-time "u" which is a reg- 
ularizing variable for only two particular bodies. It has been assumed up until now that m 0 and m, 
are the colliding bodies and the distance, r 2 * r , between the two bodies has been going towards 
zero A glance at Equations 154, will point out that for a collision between m 0 and m 2 or m, and m 2 . 
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the differential equations of motion have singular points, namely when r 0 = 0 or t x =0. There- 
fore "u" is only a satisfactory regularizing variable for two particular colliding bodies. An in- 
dependent variable "w" will be introduced that will remove all singularities and allow any number 
or combinations of collisions to occur without singularities. We define; 


dt - T dw 
t = 0 


for 


w = 0 

where 

r = (, -.-v*) 


(323) 


(324) 


Now i is defined by 291, and the function r has a determined value for each real value of time and 


0 < r < i 


(325) 


From 323, 325 it is seen thatw and t increase and decrease together. Furthermore, there 
exists a one to one correspondence between w and t . It is clear that r is positive when all the 
distances r 0 , r, , r 2 are greater than zero, and w cannot become infinite when t tends toward a 
finite value, say for t = t r From 117 and 323. 


dw r 

du ” r 


Since r , and r go to zero together, the right side of the equation remains finite for r- 0. It is also 
clear that if r - o for t ~ t , , w tends toward a finite value when u - u, ort-t,. One finds the same 
result if the distance r 0 orr, goes to zero for t = t,. The variable 'V* will be finite when t is 
finite and since 1 1 1 < |w| , from 323 and 325, the reciprocal is also true. It can therefore be said 
that 


4 

lim w - + 00 

1 im 


t= +00 

t = -00 


lim t = + 00 

1 im 


w= + 00 

W = -00 


CO 


00 
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We shall devote the next section to finding a lower limit of the radii of convergence for the co- 
ordinates of the three bodies, their mutual distances and the time expressed in power series of 
w-w. Two cases are to be investigated. 

Case 1: For w = w one of the distances r 0 , r,, r 2 is inferior to k,/2; for example 


r 


2 


< 


h 

2 


Let t, be the value of t for w = w. We shall designate by u, the value of u for t = t, or w = w. We 
can then say that the coordinates of the bodies, their mutual distances and the time will be devel- 
opable into powers of u-u, if u verifies the condition 


| u -u, | - Q 


The variables u and w are related by the equation 


du _ r 

dw ~ r 


(326) 


for 


w 


w 


Since both r and l/r are developable into a series in powers of u - u , , the right side of 326 is also 
developable into a power series if |u - u, | < Q. In order to apply the Cauchy-Picard Theorem to 
Equation 326, we must find an upper limit for \r/t\ when | u - u, | < Q. Writing down 198, 201, 184 


T o = Pi + El 

45 , 

El < 49 P i 

P, > 14k, 


we see that r 0 2 never becomes zero or negative and the real part of r 0 doesn't change sign. Since 
the real part of r 0 is positive for u = u, , it will remain positive if |u - u, | IQ. It therefore follows 
that 



< 1 
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and 



< 2 


In an analogous manner 



< 2 


Now we observe that 


1 - e' r ^ 


1 r r 2 

1 k k 2 

, 1 . 1 k i 

e k ^-l 

r 


^ ~ 2l 2 + U 3 

< -TT + + + • • < 

- * 2 V a 3 

k i 


from 293 l = 29/2 k : and e 2/29 - l < 1 / 12 , we see that 


r 

r 


1 - e' r/!t 


-'J* I 




12 


• 2-2 


1 

3k^ 


so that 


r 


r 


< 


1 

3k j 


if | u - Uj | 1 Q . 

Now we can apply the Cauchy-Picard Theorem and say that the Equation 326 


du _ r 

Hw “ T 


has a unique analytic solution (in powers of w - w) for 


| w - w] 1 3 Qkj 


(327) 


A result of this is that the coordinates of the bodies, the distances r 0 , , r and the time are 

developable into powers of w - w if w verifies 327. Since k t and Q are symmetric with respect to 
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the masses m 0 , m, and m 2 there would be no change if, in the place of r 2 = r < k,/2 one had 
r Q <k,/2 or r, <k,/2 for W = w. 

Case 2: For w = w all of the distances r 0 , r , , r are > k,/2 or from 293, r 0 , r, , r > 14k 0 . We 
note that this is condition 132 if t is the value of t for w = w . We have already found that the co- 
ordinates of the three bodies and the distances r 0 , r, , r 2 are developable into powers of t - t if t 
verifies 153. From the above reasoning 


l - 





< 2 


so that |r| <8 if 1 1 - 1 | <T' . We note from 323 that 


dt 

dw 


r 


and the Cauchy-Picard Theorem insures that this equation has a unique analytic solution in powers 
of w - w for 


|w - w| - 8 T ' 


(328) 


That is, the coordinates of the three bodies, the distances r 0 , , r 2 and the time are, in case 2, 

developable into powers of w- w if 328 is satisfied. From 153, 293 


or 


r 

8 




4M 2 

21mk r 


+ M K 


1 / 

224 y 


^ •»!« 


r« 

fXvM 7~, 

224 y 21mirj + l K l 



r 

8 


224 



k i 


From 320 


3Qk, 


15 M 3 , 

6 + TiT + 2in Gk l + 


/3k, 

TT 

9 

2n 


Gysr, 


3 


|K| k. 
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i 


Both 3Qkj andT'/8 are greater than 


( 3k , 


15 

8 




k, + 


_9 

2m 


GvMk 


| |K| 


k, + 224 




+ 3 K k, 


(329) 


where G and k, are defined by 318 and 293 and m designates the smallest of the masses m 0 , m, , m . 
An important result of Sundman’s work can finally be stated. 

The coordinates of the three bodies, their mutual distances and the time are developable in 
powers of w - w, such that these developments are convergent for 


| w - w| < fi 


We now have a "convergent strip" with a width 20 and the real w axis runs symmetrically through 
the center of this strip. Since the domain of convergence for a Taylor series is circular and we 
now have a strip of convergence, a transformation would be useful. With this in mind, we in- 
troduce a new variable r by the transformation. 


20 1 + r 

w = IT lo e F^T 


g-TT-w/20 „ j 

e 77w/2fi + | 


> 


(330) 




It can be shown that as w |r| - l and that M < i . The real values of r between -l and +1 will 
have a one to one correspondence with the real values of t between -<» and +®. We have effectively 
transformed the strip of convergence in the w plane to a unit circle in the r plane. We can now 
state Sundman's final theorem. 

Theorem 8 : If in the problem of three bodies, the constants of area are not all zero, one can 
find two constants l andfi, (the coordinates and the velocities of the bodies being given for a certain 
finite moment) such that, if one introduces in place of "t" a variable "r", the coordinates of the 
three bodies, their mutual distances and the time are developable in powers of r . These series 
solutions converge for M < i and represent the movement for all time. The equations remain 
regular for collisions between any two bodies. 

Through the introduction of the variable r, the coordinates of the three bodies, their mutual 
distances and the time is developable into powers of r if |r| < 1. In the general case considered 
by Sundman, the series are very complicated. The convergence of these series was investigated 
by Belorizky and the results are summerized below. 


INVESTIGATION OF CONVERGENCE 

After determining the constants l andfi, the new variable r is substituted for t and the regu- 
larized solution series are obtained. Belorizky uses the equilateral triangle solution to the three 
body problem in an effort to investigate the convergence of the Sundman solutions. He takes the 
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mutual distances between the three bodies as one astronomical unit, the sum of the masses is the 
unit mass and the Gaussian constant is set equal to one. Using these simplifications r is a constant 
and 


r 



consequently t - Tw and 


t 


2ar 1 + r 

7T 1 - T 


1 + T 

A log x ~ 


and 


t 


2A r + -3 



+ 


r 2n + l 

+ 2iT+l + 


Furthermore 


X - cos t 


t 2 t 4 

1 - 21 + 4T - •" 


2A 2 t 2 - ( 


3*’ 


A 4 ) 


(S**- 


A 4 + 


*»•)■ 


sin t - t 


5! 


= 2Ar + (^ - J A 3 ) r 3 + (^ - J A 3 + ^ A 5 ) r 3 + 


We consider m 2 to be the larger of the three masses, and m, - m 0 . If 

1) m i = 200 

2 ) "h - K) 

3) ">i = 3 
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the corresponding values of f) will be; 

1) n, < 9 x io" 8 

2) n 2 < 4 x io - 6 

3) n 3 < 1(T 5 

so that A has the values 

1) Aj < 4 X 10~ 8 

2) A 2 < 2 x IO' 6 

3) A 3 < 4 x 10' 6 

Let "h" be the error introduced by stopping the "t" series at 2n - 1 terms so that 

( t 3 t 5 r 2 "" 1 

t - h - 2A^r + -y + T + • • • + i 

and x, y can be expressed as 



1 L 

/ t 3 r 5 

T 2n-1\ 

i 2 1 r 

( t3 

X = 1 - 

• 2T [2 

A V + T + T + '•* 

+ 2n- 1/ 


k \ T + ~S + • • 

1 

' T 3 

r 5 r 2n-l 

\ i r 

/ T 3 

r 2n ~l \ 

< 

CN 

II 

>> 

v T + 3 

+ 5 + ‘ ‘ + 2n - 1 

)"3T [ 2 

A ( t+ ~3 + ■ 

* * + 2n - 1 / 


If we wish to have the coordinates x, y for the epoch t = 1 with an accuracy of only one decimal 
place, the number of terms necessary (n) can be computed by setting 


1 - 2A 


(-3 


+ 5" + • ' • + 2n 




101 



hence 


L 


( t 3 T 2n-1\ 

2A\r + T 2^l] 


> 0.9 


and the sum 


T 3 T S T 2 "' 1 9 

S r " T + T + T + '" + 2 FT > 20S 


The variable r is determined by the relation r - e 1/A - l/ e 1/A + 1 and for 


A<4 x 10~ 8 on finds S Tj >10 7 
A < 2 x 10' 6 on finds S T > 2 x 10 s 
A < 4 x 10' 6 on finds S T >10 s 

T 3 

We shall now compare the sum S T with the divergent series 


1 


+ 



1 

7 + 


1 

+ 2n- 1 


+ ••• 


Since r < 1 one has r 2 "* 1 / 2p- 1 < l/2p- 1 and if 


S = £ 2p-l 

1 

Belorizky shows that p > e 2s_2 - 4 and since S T <s, p> e T . Now for 

S T = 10 7 one finds p, > 10 8X1 ° 6 

S = 2 x 10 s one finds P , > 10 ,7X1 ° 4 

t 2 1 

S = 10 s one finds p, > 10 8X10< 

T 3 3 

We have designated by "n" the number of terms which is necessary in the series 

t 3 r 2n_1 _ J_ 

T +^-+---+ 2n - i + "- - 2A 


102 


/ 


in order to have S T > 9/20A. This number n is greater than the number "p" of the terms which it is 
necessary to take in the series 


. 1 1 1 

1 + 3 + 5 + 7 + 


in order to have the same sum. From above we can see thatn > p. 

We wish to calculate x, y in our particular case with only a single decimal accuracy by Sundman's 
method. For the epoch t = 1, in the three different cases where M/m = 200, M/m = 10, M/m = 3, the 
number of terms is respectively greater than IO 8 * 106 , 10 17X1 ° 4 , 10 8X104 . 

The extreme slowness of convergence in this particular case is apparent. The Sundman solution 
is, in a practical sense, quite useless for ephemeris computation. The motion of a body near col- 
lision, can however be represented by the Sundman series solution. The Sundman exposition, al- 
though not generally useful for a solution to the three body problem, has allowed a rigorous in- 
vestigation of the motion near collision. A few of the more important results are given below. 


SUMMARY 


1. In a double collision, the distance to the non-participating third body remains bounded. 

2. In order for a triple collision to occur, the area integral of the system must equal zero. 

3. The velocity and acceleration of the colliding body tends toward infinity at the point of 
double collision. 


4. The unit vector along the radius joining the two colliding bodies tends to a limit. 

5. Singularities in the equations of motion can be removed at one of the bodies by the intro- 
duction of the independent variable "u". 

6. An independent variable "w” was found such that all the singularities of the motion were 
removed, allowing any number or combinations of double collisions to occur. 

7. The power series solutions were found to converge, although extremely slowly. 

Although this paper, up until this point, has considered only the regularization introduced by 
Sundman, there are several other known transformations that regularize the equations of motion. 
One of the best known methods for the restricted three body problem is that of Levi-Civita. 

We consider two primary bodies in the cartesian plane having the coordinates (/x, 0), (m - 1, 0). 
The third body's position is denoted by (x, y). A transformation is made from the x, y plane to a 
p, q plane by the following equations 


2 2 
x “ M - P “ q' 


y - 2pq 

dt = 4(p 2 + q 2 ) dr 
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Using these transformations, one can transform the equations of motion in the (x, y) plane 


d 2 X 

dt 2 


dy 
2 dx 


n 


X 


d 2 y 

dt 2 


+ 


dx 


2 at 


n 


y 


n 


\ [(1 -M) r 2 +Mr 2 2 ] 


+ 


1-M 


+ 


JL 

r 2 


rj = )/(x-m ) 2 + y 2 


r 2 = /(x + 1 -M) 2 + y 2 


to the corresponding equations of motion in the (p, q) plane. These are; 


^7 * 8(p j + q 2 ) af = [4(0- 2 c)(p J + q 1 )] 


(£)'*(£)’ = 8(0-^c)(p. t q>) 


Subscripts in the right side of the equations denote differentiation with respect to that variable 
and C is a constant of integration. 

If we now consider a collision at (^, 0), then as r, - 0 the term (l-rf/r, or (1 - M )/ P 2 + q 2 in 
fi which becomes infinite, appears multiplied by a factor p 2 + q 2 . At p - q - 0, the last equation 
shows that the square of the velocity is proportional to 8( 1 - M ) . The curves of motion are analytic 
curves without singularity near the origin. Thus for a single collision at (m, 0), the equations are 
regular. 

Recently, a new method has been introduced into the three body problem. Originating in atomic 
physics "spinors" were used to describe the process of spinning of an elementary particle. In 
celestial mechanics, spinors are no longer used to describe spin but are used instead as a mathemat- 
ical aid in simplifying the equations of motion for the three body problem. A spinor can be thought 
of as a vector in the complex plane connecting a complex number *, = u, + iv, to another complex 
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number z 2 = u 2 + iv 2 . The spinor is then a four space vectors of the form: 



It is these vectorial analogies of spinors which are used to regularize the differential equations of 
motion. Transformations of the Levi-Civita type cannot immediately be extended into three di- 
mentions. By using spinor notation, a position vector r and its 3 scalar components (x, y , z) can 
be expressed in terms of z , , z 2 and their complex conjugates. 

For each position vector, there corresponds differential equations for a spinor (z„ z 2 ) which 
are regular, when a certain pseudo-time is chosen to be the independent variable. 
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